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Abstract 

I derive a Wick ordered continuous renormalization group equation for 
fcrmion systems and show that a determinant bound applies directly 
to this equation. This removes factorials in the recursive equation for 
the Green functions, and thus improves the combinatorial behaviour. 
The form of the equation is also ideal for the investigation of many- 
fcrmion systems, where the propagator is singular on a surface. For 
these systems, I define a criterion for Fermi liquid behaviour which 
applies at positive temperatures. As a first step towards establishing 
such behaviour in d > 2, I prove basic regularity properties of the 
interacting Fermi surface to all orders in a skeleton expansion. The 
proof is a considerable simplification of previous ones. 







1 Introduction 



In this paper, I begin a study of fermionic quantum field theory by a contin- 
uous Wick ordered renormalization group equation (RGE). As an example, 
I take the standard many-fermion system of solid state quantum field the- 
ory, but the method applies to general fermionic models with short-range 
interactions. I show that a determinant of propagators appears in the RGE 
and I use a determinant bound to prove that a factorial which would appear 
in bosonic theories is removed from the recursion for the fermionic Green 
functions. This may lead to convergence of perturbation theory in the ab- 
sence of relevant couplings, but I do not address the convergence problem, 
which is related to the solution of a particular combinatorial recursion, in 
this paper. A short account of this work has appeared in J13|], 

Continuous RGEs were invented by Wegner |] and Wilson @. Polchin- 
ski ||] found a beautiful way to use them for a proof of perturbative renor- 
malizability of c/> 4 theory. His method was simplified in ||, and extended 
to composite operator renormalization and to gauge theories by Keller and 
Kopper ||. Keller j7| also proved local Borel summability. While equiva- 
lent to the Gallavotti-Nicolo ^, [n| method, the continuous RGE is much 
simpler technically. An application of continuous RG methods to nonpertur- 



bative bosonic problems [11, [T^] requires many new ideas and a combination 
with cluster expansion techniques, to control the combinatorics. It is one of 
the points of this paper that for fermions, the straightforward adaption of 
the method yields a determinant bound which improves the combinatorics 
of fermionic theories as compared to bosonic ones, and may lead to nonper- 
turbative bounds. The determinant structure is not visible in the form of 
the flow equation used in || ^, because the flow equation in that form 
is a one- loop equation which has too little structure. A key ingredient for 
the present analysis is Wick ordering, which was first used in the context 
of continuous RGEs for scalar field theories by Wieczerkowski Q. I show 
that for fermions, the Wick ordered RGE contains a determinant of propa- 
gators to which a Gram inequality applies directly. A closer look at the way 
the Feynman graph expansion is generated by the RGE shows that the sign 
cancellations bring the combinatorial factors for fermions nearer to that of 
a planar field theory. This reduction does, however, not lead to a planar 
field theory in the strict sense because of certain binomial factors in the 
recursion. 

The direct application of the determinant bound shown here requires 
the interaction between the fermions to be short-range. This prevents a 
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straightforward application to systems with abelian gauge fields by simply 
integrating over the gauge fields. One model to which the method applies 
directly is the Gross-Neveu model (which has been constructed rigorously 



4], 15]). I show here only the most basic power counting bounds by leaving 



out all relevant and marginal couplings, but it is possible to take them into 
account by renormalization. 

A class of physically realistic models with a short-range interaction is 
that of nonrelativistic many-fermion models. In these models, there is a 
significant complication of the analysis because the singularity of the fermion 
propagator in momentum space is not at a point, but instead on the Fermi 
surface, which is a (d — l)-dimensional subset of momentum space. Only in 
one dimension, the singularity is pointlike - the 'surface' becomes a point. 

The interest in these models has resurged recently because of the dis- 
covery of high-temperature superconductivity. Before that it was taken for 
granted that Fermi liquid (FL) behaviour holds in all dimensions d > 2, 
and Luttinger liquid behaviour in one dimension (the latter has been proven 
[fi6| , [[?], [l8|]). At certain doping values, however, strong deviations from FL 
behaviour are seen in the high-T c materials. The discussion following these 
discoveries revealed that the former arguments for FL behaviour contained 



logical gaps. Like [19, 2C, 21, |23|, 24, 25], the present work is not aimed 
at an understanding of these deviations, but at the more modest goal of 
determining first under which conditions FL behaviour occurs. 

Before doing so, it is necessary to give a definition of what would con- 
stitute FL behaviour. At zero temperature, the noninteracting Fermi gas 
has a discontinuity in the occupation number density; this is also a property 
one would require of a zero temperature FL. This discontinuity is absent 
for Luttinger liquids. It is, however, not sufficient for FL behaviour because 
there is a one-dimensional model which has both such a discontinuity and 
some Luttinger liquid features in the spectral density ]2q| . Moreover, in the 
standard models of many-fermion systems such a step never occurs because 
superconductivity sets in below a critical temperature, and it smoothes out 
the step in the zero-temperature Fermi distribution. It is thus desirable to 
give a definition of FL behaviour at temperatures above the critical tem- 
perature for superconductivity. This is not at all straightforward because 
there is no clean characterization of FL behaviour at a fixed temperature. I 
propose to look at a whole range of temperatures and values of the coupling 
constant to bring out the characteristic features of a FL. 

I define an equilibrium Fermi liquid as a system in which the perturba- 
tion expansion in the coupling constant A converges for the skeleton Green 
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functions in the region |A| log/3 small enough (here (3 is the inverse temper- 
ature) and where the self-energy fulfils certain regularity conditions. The 
skeleton Green functions are defined in detail in Section ^; in them, self- 
energy insertions are left out, so that the Fermi surface stays fixed. I discuss 
in Section |7] how they are related to the exact Green functions. The logarith- 
mic dependence of the radius of convergence on (3 comes from the Cooper 
instability; the difference between Fermi liquids and Luttinger liquids is in 
the regularity properties of the self-energy. This is discussed in detail in 



Section 2.6 



The goal is to show that the standard many-fermion systems are Fermi 
liquids in that sense. A proof of this requires a combination of the regularity 
techniques of [21, 23] for renormalization with the sector technique of 



20| in the determinant bound (see Section 2J3 for further discussion). I do 
not give a complete proof in this paper but only a part of it, by showing a 
determinant bound and some of the required regularity properties of the self- 
energy in perturbation theory. The hope is that the determinant bound will 
lead to convergence, so that the method developed here, which is somewhat 
simpler than e.g. the one in [f20| , will work nonperturbatively. A different 
representation for fermionic Green functions that provides a simplification 



and leads to nonperturbative bounds is given in [27]. 

In Section I review the Grassmann integral for many-fermion systems 
briefly, to give a self-contained motivation for the study of such systems, 
and to fix notation. The Fermi liquid criterion is formulated in Section 



2.6| . Section |3j contains the general renormalization group equation and the 
determinant formula Eq. (|53|). Section |] contains the determinant bound 
and an application to systems where the propagator has point singularities. 
In Section [|, I show the existence of the thermodynamic limit for the many- 
fermion system in perturbation theory. In Section ^, I prove bounds on 
the skeleton self-energy that are needed to renormalize the full theory in 
perturbation theory. Again, the RGE in the form of ||, ||, |^] would not be 
very convenient for this because it is a one-loop equation, whereas the crucial 
effects for regularity all start at two loops. They can be seen in a simple way 
in the Wick ordered RGE. Details about Wick ordering and the derivation 
of the determinant formula Eq. (|53|) are deferred to the Appendix. 



Acknowledgement: I thank Volker Bach, Walter Metzner, Erhard Seiler, 
and Christian Wieczerkowski for discussions. I also thank Christian Lang 
for his hospitality at a very pleasant visit to the University of Graz, where 
this work was started. 
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2 Many-fermion systems 



The model is defined on a spatial lattice with spacing e. Continuum models 
are obtained in the limit e — ► 0; lattice models, such as the Hubbard model, 
are obtained by fixing e. 

Let d > 2 be the spatial dimension, e > 0, L 6 R be such that | £ N, 
and let G be any lattice of maximal rank in R rf , e.g., G = Z d . Let L be the 
torus L = eG/LG. The number of points of this lattice is |L| = (j) d - Let 
/ L dx F(x) = £ d £ x6L F(x) and «5 L (x,x') = e~ d 5 x y. 

Let be the Fock space generated by the spin one half fermion op- 



erators satisfying the canonical anticommutation relations [28], i.e. for all 
x, x' E L, 

c a( x )c+,(x') + C a'( X ') c a( x ) = <W ^l(x,X ). (1) 

Here a E {—1,1} is the spin of the fermion in units of |. The free part of 
the Hamiltonian i3x,(c, c + ) = Hq + XV is 

H = - £ /dx /dyT(x,y)c+(x)c a (y). (2) 

For a one-band model on a lattice with fixed spacing e, T(x, y) = i x -y = 
iy-x describes hopping from a site y to another site x with an amplitude 

t x — y — ^y— x- 

The interaction is multiplied by a small coupling constant A; I assume it 
to be a normal ordered density-density interaction 

V(c,c + ) = - J dx J dy Yl «( x -y)4(x)c+(y) CQ (x) C(T (y) (3) 

L L a,ae{-l,l} 

In other words, it is a special type of a four-fermion interaction. 

For instance, the simplest Hubbard model is given by A = §■ where U is 
the usual Hubbard-f/, and by v(pt — y) = <5l(x, y), and the hopping term is 
t x _ y = t if |x — y | = 1 and zero otherwise, where t is the hopping parameter. 

At temperature T and chemical potential fi, the grand canonical parti- 
tion function is given by Zl = tr^ L (q-P( h ^-v n i-)} with Nj^ = J h dx n(x), 
and (3 = j^ji- Observables are given by expectation values of functions, 
mainly polynomials, of the c and c + , 

(0) h = ±- tr (e-^-^)O(c, c+)) (4) 
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A basic question is whether the expected values of observables have a finite 
thermodynamic limit and whether an expansion in A can be used to get their 
behaviour at small or zero temperature T. For instance, one would like to 
expand the two-point function (c + (x)c(v))l = Yl^Lo -^^2 r ( x ' y)- ^ ^ s by 
now well-known that the result of a naive expansion e~ xv in powers of A is 



that at T = 0, lim^oo = oo for all r > 3 (see, e.g., [|19| |2lJ ) . At positive 
temperature T, this unrenormalized expansion converges for |A| < const T d ; 
see Section [|. To get a better T-dependence of the radius of convergence, one 
has to renormalize. Because of the BCS instability, the best one can hope 
for in general is a bound |A|log^ < const for the region of convergence. 
This is part of the Fermi liquid criterion formulated below. 

2.1 Grassmann integral representation 

The standard Grassmann integral representation is obtained by applying 
the Lie product formula 



o 



n T — >oo \ / 



to the trace for Z A and (O). The spacing in imaginary-time direction is 
e t = The limit exists in operator norm because on the finite lattice L, 
all operators are just finite-dimensional matrices. 

Inserting the orthonormal basis of J\, between the factors in Eq. (||) and 
rearranging, I get Z L = lim^^ Z LiTIt , where Z Lj „ t is given by a finite- 
dimensional Grassmann integral, as follows. Let n T be even and T = {t = 
ne T : n G Z, -^f < n < ^}, let A = T x L, and A be the Grassmann 
algebra generated by tp a (x), 4> a (x), with a € {1,-1} and x = (r, x) € 
A. Fix some ordering on A and denote the usual Grassmann measure p9| , 
H by D A ^D A ^ = \[ X (T d^ a {x) d$ a (x). Then Eq. (|) implies Z Li „ T = 

N A j D A ^D A $ e ~ s ^'^ where A/a is a normalization factor that depends 
on e,L, and n T , and where 

S A (^) = j T dT {Y,J ^(r,x)a T ^(r,x) - H A (t/j(r),^(r))\ (6) 

Here I have used the notations ^(r)(x) = ip(r, x), J" T dTF(r) = e T J^tgt F( t )i 
and d T Tp(r) = e~ (ip a (r+£ T ) — t/> ct (t)), and the sum over r runs over T, with 
antiperiodic boundary conditions |pH|. For n T < oo and L < oo, this is a 
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finite-dimensional Grassmann integral. The limit n T — > oo, and afterwards 
L — > oo, will be taken only for the effective action. No infinite-dimensional 
Grassmann integration will be required. 

To do the Fourier transformation, it will be convenient to deal with 
periodic functions defined on an interval of double length in r, and to impose 
the antiperiodicity as an antisymmetry condition: let T2 = e T Z/2/3Z, in 
other words, T2 = {r € e r Z : — (3 < r < /3} with periodic boundary 
conditions. Thus the fields ip and ip are periodic with respect to translations 
of r by 2/3, and antiperiodicity with respect to translations by (3 is imposed 
by setting 

V(r + /?,x) = -V(r,x) (7) 

for all x £ A. With the further notation J* A dx F(x) = Jjdr J* l dx F(t, x), 
and 5\ (x, x') = e T ~ 1 e~ d 5 TtT i5 XjX i , the action is S(ip, V>) = 82(1^, fy+XS^tp, tp) 
where 




with 

a(x, a, x', a') = S aa > ((d T + n)S A (x, x') - T(x, x')5t 2 (t, t')) (9) 

and 

S^^) = Y, j dx j dx' ^(x)^ a (x)v(r,x,T f y)$ a ,{x')^,{x') (10) 

<T ' <T ' A A 

with v(t, x, r', x') = <5t 2 (r, r')v(x — x'). For the present work, the interaction 
does not have to be instantaneous. Retardation effects, like from phonons, 
are allowed. That is, v(t, x, t', x') may have a dependence on r and need 
not be local in r. 

The operator a appearing in S2 is invertible because the antiperiodicity 
condition removes the zero modes of the discretized time derivative. In 
other words, the Matsubara frequencies for fermions are nonzero at positive 
temperature (this will become explicit in the next section). 
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2.2 The propagator in Fourier space 

Fourier transformation with the antiperiodicity conditions Eq. (|7|) is de- 
scribed in Appendix [A|. The Fourier transforms of ip and $ are 

Mp) = J dxe-^Mx) (11) 

A 

where, for p = (cj, p) and x = (r, x), px = lot + px. If L* is the dual lattice 
to L, the momentum p is in A* = M„ r x L*, where 

7T Th Tl 

M nr = K = ^(2n + l):neZ, -y<i<y} (12) 

is the set of Matsubara frequencies u n . With the notation J* A „ dp F(p) = 
^E^ eMnr k* d P F (v,p), where / L , dp = L~ d £ peL ,, the inverse Fourier 
transform is ip a {x) = J A » dx e ipx ip a (p). The Fourier transform of the hop- 
ping term is T(p, q) = 5^* (p + q, O)T(q), with T(q) = f dz e iqz i z . Denoting 

L 

E(p) = f (p) - /x, (13) 
where /U is the chemical potential, 

oij = t^— (e iE ^ - 1) , (14) 

and (5a*(p + p' ; 0) = <5l(p + p',0) er^^-wy, the Fourier transform of the 
operator a in the quadratic part of the action is 

a(p,a,p , ,a') = S A ( P + p',0)S aal (Q - E(p'j) (15) 



In other words, the matrix with entries a(p,cr,—p',cr') is diagonal, and for 
temperature T = i > 0, all diagonal entries are nonzero because 



Re uj\ 



sin(e r a;) 



e T uj 

,-i 



a w < 16 > 



Thus the inverse of a, the propagator c = a , exists; it has the Fourier 
transform 

c(p, a, p, a') = (5 A * (p + p', 0)<W — — (17) 

iu' — E(p') 
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In the formal continuum limit e T — > 0, Q — > u, so one gets the usual formula 
(ioj — E(p))~ l . The partition function of the system of independent fermions 
(A = 0) is 

D A iP D A ij} e^' A ^ = det A = JJ (m>) - E(p)) (18) 



which is nonzero by Eq. (16). 



2.3 The class of models 

Denote the dual to G by B, the first Brillouin zone of the infinite lattice. 
For instance, for G = eZ rf , B = W 1 /^-!/ 1 . The assumptions for the class of 
models are: there is ko > 2 such that the dispersion relation E £ C k °(B,M,), 
and for all p£B, E(— p) = -E(p) holds. The interaction v is a C k ° function 
from M x B to 1R, all its derivatives up to order ko are bounded functions on 
SxM, v(— po, p) = u(po) p)> and the limit po — > oo of exists and is C fc ° in p. 
There is go > such that for all p on the Fermi surface S = {p : E(p) = 0}, 
|V-E(p)| > go holds. The Fermi surface is a subset of an e-independent 
bounded region of momentum space (hence compact), it is strictly convex 
and has positive curvature everywhere. In particular, there is Vi > such 
that for all L and e 

J dkl(\E(k)\ < 2) < Vi. (19) 

The constant 

£max = SU P |£(p)| (20) 

pee 

is independent of e. 

Under these hypotheses, there is eo > and a C 2 -diffeomorphism % from 
(-2e ,2e ) x 5 d_1 to an open neighbourhood of the Fermi surface S in B, 
(p, 0) i— > w(p, 9), such that 

E(w(p,6))=p and \d p *(p, 0)\ < — (21) 

9o 



(sec |2T]], Lemma 2.1, and [22|, Section 2.2; eo was called ro there). Let 



J(p,9) = det nr' (p, 9) and denote 

J = sup |J(p,0)| and J x = sup / d9\J(p,9)\. (22) 

IpI< £ loKen J 



p\<<-{> |p|< e o 
eesd- 1 S d ~ 
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I assume that eo < 1 and (for convenience in stating some bounds) that 
/3eo = -j^f > 6. With the units chosen in a natural way, i.e., with typical 
bandwidths of electron volts, this corresponds to temperatures T up to 1000 
Kelvin if eo is of order one, which seems a sufficient temperature range to 
study conduction in crystals. Note, however, that eo depends on the Fermi 
surface and thus on the filling factor. A typical example is the discretized 
Laplacian 

1 d 

E(k) = S-2^(l-cos(efc 1/ ))-Ai. (23) 

£ u=l 

For e — > and B = — , E(k) — > k 2 /2m — the Jellium dispersion relation, 
which satisfies the above hypotheses if large |k| are cut off. For e = 1 
and B = t, one gets the tight-binding dispersion relation with hopping 
parameter t/2, which satisfies the above hypotheses if n ^ td (half-filling) . 
In the limit \i — > td, eo — > in the Hubbard model. This implies that to 
have bounds uniform in the filling, one has to stay away from half-filling. 
The energy eo sets the scale where the low-energy behaviour sets in. The 
effective four-point interaction at that scale can differ substantially from the 
original interaction. For a discussion, see Section R and [24]. 



2.4 Nambu formalism 

It will be useful for deriving the component form of the RGE to rename 
the Grassmann variables such that the distinction between ip and V> is in 
another index. This is a variant of the usual 'Nambu formalism'; see, e.g., 
g§. Let 

T = Ax {-1,1} x {1,2}, (24) 

and denote X = (x,a,i) E T. For x E A and a E {—1,1}, the fields are 
defined as ip(x,a,l) = tp a (x) and tf)(x,a,2) = ipa(x). The antiperiodic- 
ity condition reads ip(x + (3e T ,a,i) = —ip(x,cr,i) with e T the unit vector 
in r-direction. The Grassmann algebra generated by the (ip(X))xer is de- 
noted by -4.r[V']- Given another set of Grassmann variables (r](X))xeVi 
the Grassmann algebra generated by the ifi and r\ is denoted by ^IrfV') 7 ?]- 
Furthermore, denote J r dX F(X) = Yli=i Ylae{-1 1} Ja ^ x cr, i) and 
Sr((x, a, i), (a/, a', i')) = 8n> 5 acr i 6\(x,x'), and define a bilinear form on 
A T [ip,v] by 

(^ n) r = f dX ij(X) n(X) = - ( V , . (25) 
r 
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Then S2 = \ (ip, A ?/>) r where, for X = (x, a, i) and X' = (a/, a 
A(X,X') 



' -> i') 



( if i = i' 

a(x, a, x' , a') if i = 1 and i' = 2 (26) 
— a(x', cr', x, a) if i = 2 and i' = 1 



with a given by Eq. (g). In other words, when written as a matrix in the 
index i, A takes the form 

a 
-a T 



(27) 



with (a T )(2;, cr, x', a') = a(x', a', x, a) denoting the transpose of a. Since a is 
invertible, A is invertible as well. 

With this, Z L ,n T = JVAdeta where Z A = J dfi C We~ XS4W , where 
C = A -1 , and d/ic 1 is the linear functional ('Grassmann Gaussian measure') 
defined by dycW = (deta) -1 D T tp e^'^ . The constant N A det a drops 
out of all correlation functions and can therefore be omitted. The 'measure' 
due is normalized, J dyci'P) = 1; an d its characteristic function is 

^(^eW'^e^^*. (28) 

All moments of dfic can be obtained by differentiating Eq. (|28|) with respect 
to rj and setting r\ = 0; see also the next subsection. 

2.5 The connected Green functions 

In the correspondence between the system, as defined by the Hamiltonian 
Hi, and to the Grassmann integral, I have so far only discussed the 
partition function itself. In the path integral representation of Eq. (Q), with 
a polynomial observable 0(c,c + ), one simply gets a factor O(ip(0),ip(0)) 
in the Grassmann integral. The m-point Green functions of the system 
determined by C and V are 

I m \ ~\ r m 

/ I] j,(X k ) ) = jrj ducm e- xv ^H^(X k ). (29) 

\ k == 1 / k == 1 

They determine the expected values of all polynomials by linearity. 

On a finite lattice, the limit n T — > oo of Z\ exists by the Lie product 
formula Eq. (^j, and for A small enough (depending on L, e, 0, and y), it is 
nonzero since the trace of the matrix Q-Pi H ^-^ N ^) over the finite-dimensional 
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space T\, is a continuous function of A, which is nonzero at A = by Eq. 
( |l6|) and Eq. (|l~8| ) . A similar argument applies to the numerator of Eq. (f|) . 
Thus the limits of numerator and denominator in Eq. as n T — > oo exist 
separately. Therefore one can take this limit in numerator and denomina- 
tor through the same sequence, i.e., take n T to be the same in numerator 
and denominator. It follows that with the special choice O(ip(0),ip(0)), all 
expectation values in the Hamiltonian picture can be expressed as the limit 
n T — > oo of Eq. (p9[), with a special choice of the polynomial in the fields. 
Thus the correlation functions given by Eq. (|2^) include as a special case the 
expectation values of polynomials in the creation and annihilation operators. 

Let (r](X))xer be a family of Grassmann generators. The partition func- 
tion with source terms is Zr(r)) = J dfj,c(ip) e~ xv ^ + (*?'^) r . Let 



-1 c -d d 
- £ d V (X) 



then 



n ^ 



\k=l 



Z T (0) 



n 



(30) 



r?=0 



Thus, if one knows Zi>(rj) one can derive all correlation functions. 

It is convenient to study the connected correlation functions, defined as 



\fc=i 



n 



log Z v {r\) 



(31) 



rj=0 



instead. Since Z-p(r]) is the exponential of log Zr(rj), one can reconstruct all 
correlation functions from the connected ones. 

It is even more convenient to transform the sources rj, to get the ampu- 
tated connected Green functions. They are generated by 



Q eS (x) = log / dfidij) e - AV ^+*) 



(32) 



A shift in the measure shows that G e ff(x) = \ (x 5 C~ 1 x) r + log Z? (C _1 x) • 
so that the study of Q c s is equivalent to that of log Zy ■ 

The self energy S(p) is defined as the one-particle irreducible part of the 
two-point function. In terms of the connected amputated two-point Green 
function G*2, which is the coefficient of the quadratic part (in if)) of the 
effective action Q e s(ip), it is E(p) = G2(p)(l — CG^p)) -1 . 
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2.6 Criteria for Fermi liquid behaviour 

In the following I give a definition of Fermi liquid behaviour which is linked 
to the question of convergence of the expansion in the coupling constant 
A, and I discuss in some detail the physical motivation for this definition, 
the results that have been proven in this direction, and its relation to other 
notions of FL behaviour. 

In most many-fermion models, one cannot expect the expansion in A to 
converge uniformly in the temperature, not even after renormalization. In 
particular, the Cooper instability produces a superconducting ground state, 
and thus a nonanalyticity in A, if the temperature is low enough. This hap- 



pens even if the initial interaction is repulsive 34]. Nesting instabilities 
can produce other types of symmetry breaking, such as antiferromagnetic 
ordering, which may compete or coexist with superconductivity, but the 
conditions I posed, in particular the curvature of the Fermi surface, remove 
these instabilities at low temperatures (which temperatures are 'low' de- 
pends on the scale 6q). 

Let the skeleton Green functions be defined as the connected amputated 
m-point correlation functions where self-energy insertions are left out. These 
functions are the solution of a natural truncation of the renormalization 
group equation; they are defined precisely in Section ||. In particular, the 
skeleton selfenergy S s k is the second Legendre transform of the two-point 
function. 

Definition 1 The d- dimensional many-fermion system with dispersion re- 
lation E and interaction V shows (equilibrium) Fermi liquid behaviour if the 
thermodynamic limit of the Green functions exists for |A| < Ao(/3), and if 
there are constants Mq, Mi, Mi > (independent of (3 and X), such that the 
following holds. The perturbation expansion for the skeleton Green func- 
tions converges for all (A, (3) with | A | log /5 < Mq, and for all (A,/3) with 
| A | log /3 < the skeleton self-energy X s k :lxB^C satisfies the regu- 
larity conditions 

1. E s k is twice differentiable in p and 

max\\8 a E at \\ ao <M 1 (33) 

\a\=2 

2. the restriction to the Fermi surface ^sk\{o}xS € C k °(S,M.), and 



max ll^Sskl^ < M 2 . (34) 

\a\=k 
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Here k$ > d is the degree of differentiability of the dispersion relation 
E (given in Section \2l\). 



Nothing is special about the factor ^ in the condition |A| log/5 < One 
could instead also have taken any fixed compact subset of {z : \z\ < Mq}. 
The derivatives mean, when taken in po, a difference ^(E s k(po + ^jfjP) ~~ 
^sk(pO)P))- The maximum runs over all multiindices a € Nq +1 . 

This definition only concerns equilibrium properties of Fermi liquid be- 
haviour; it does not touch phenomena like zero sound, which require an 
analysis of the response to perturbations that depend on real time. It is 
natural in that it defines a Fermi liquid above the critical temperature for 
superconductance: at a given A, the value of T for which the convergence 
breaks down is T c oc e _M °/' A ', which is the usual BCS formula. Convergence 
of perturbation theory above T c implies that the usual Fermi liquid formulas 
are valid there. Convergence is stated only for skeleton quantities because 
that is all one can show. This convergence and the regularity properties 
of the self-energy, imply that the exact Green functions (no restriction to 
skeletons) are continuous in A, that the exact selfenergy £ is C 2 in A and 
p, and that £ obeys a bound similar to Eq. ((3*|). The Green functions are 
not analytic in A because otherwise already the unrenormalized expansion, 
which diverges termwise, would converge. The regularity properties (1) and 
(2) ensure that the exact Green functions can be reconstructed from the 
skeleton Green functions by renormalization. The usual skeleton expansion 
argument fl35fl , where finiteness only of the skeleton self-energy, but not of its 
derivatives, is shown, is insufficient to do that; one has to prove regularity 



properties (1) and (2). This was discussed in detail in [22|, see also Section 
0- 

The condition ko > d is necessary to make the regularized propagator 
summable in position space. It is required in the proof of Lemma || (and in 
the proofs in [20|, only that there the dispersion relation was taken C°°). In 
the absence of level crossing, the free dispersion relation E(k) is usually even 
real analytic in k. However, when reconstructing the exact Green functions 
from the skeleton Green functions, one needs regularity of the dispersion 
relation of the interacting system, and thus regularity property (2), which 
is rather hard to verify even in perturbation theory. Thus it is desirable to 
take the smallest possible k$ > d. 

Because £ obeys a bound similar to Eq. (|33|), one can do the usual 
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first-order Taylor expansion in the momenta to get 

E(p) = p (3 £)(0, P(p)) + (P - P(p)) • VS(0, P(p)) + E(p) (35) 

from which one obtains a finite wave function renormalization 

Z{p) = l + i(«9 £)(0,P(p)) (36) 

and a finite correction to the Fermi velocity, and the Taylor remainder van- 
ishes quadratically in the distance of the momentum (pq, p) to its projection 
(0,P(p)) to the Fermi surface. 

This property distinguishes Fermi liquids from other possible states of 
the many-fermion system, such as Luttinger liquids: In one dimension (where 



'Luttinger liquid behaviour' has been proven |16|, [T^, |T^] ) , the second deriva- 
tive of even the second order skeleton selfenergy grows like (3 for large (3 and 
thus violates the condition that the second derivative should be bounded 
independently of (3 for |A| log f3 < ^rp. Note that this distinction can only be 
made if [3 is allowed to vary; at fixed f3, the requirement that something is 
bounded independently of (3 is trivial. This is the reason why a whole range 
of values of (3 and A is included in Definition [j]. 

A full proof that the models obeying the hypotheses stated in [2J^ are 
Fermi liquids in the sense of Definition |l] is not within the scope of this 
paper, but several ingredients for such a proof are already in place. 

The convergence of the skeleton expansion follows for d = 2 spatial 
dimensions from a modification of the method in |20| |. The required modifi- 
cation is to put in the four-point functions. This is not difficult because at 
positive temperature, these functions have no singularities, but are bounded 
(in momentum space) by const log (3. If the recursion Eq. (|85| ) has an expo- 
nentially bounded solution, the method developed here can also be extended 
to prove this convergence. 

Regularity property (1) is proven in Section ^. Regularity property (2) 
was proven in two dimensions in perturbation theory for ko = 2 + h, h < ^ 
pl| . 22, The methods developed in Section |6] provide a simplification of 



most of these proofs, and can also be extended to give a full proof of (2). 

There is one case in which (2) becomes trivial: for the Jellium dispersion 
relation k 2 /2m— fi, is a constant, hence C°°. In that case, it is possible to 
establish Fermi liquid behaviour by a combination of the determinant bound 
of Section 0, the sector decomposition of p0| , and the overlapping loop 
method of fl2lf, without the major complications of the regularity proofs. 
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A natural question is if there are criteria independent of temperature 
that can be applied also in the zero temperature limit for 'Fermi liquid 
behaviour'. For the class of models specified in Section the simplest 
criterion is that the non-ladder skeleton Green functions are analytic in the 
coupling constant A, that the non-ladder skeleton selfenergy is C 1 , and 
that £W[5 is C k °, with bounds uniform in /3. The non-ladder skeleton 
Green functions are obtained by removing all ladder contributions (defined 
in Section ^) to the Green functions. The regularity implies that the Fermi 
velocity and the wave function renormalization are finite uniformly in /3, and 
even at zero temperature, for d > 2 (which is the usual criterion for Fermi 
liquids), whereas they still diverge in one dimension as (3 — > oo [|l^, 17, 18 1. 

The Fermi liquid criterion given in Definition [l] is more natural than the 
one using the non-ladder Green functions: if the regularity conditions (1) 
and (2) of Definition p] hold, the transition from the exact Green functions 
to the skeleton Green functions is a matter of convenience, but the replace- 
ment of the skeleton Green functions by the non-ladder skeleton Green func- 
tions changes the model drastically because it removes superconductivity. 
Evidently, a definition referring to a modified model is not as natural. Mo- 
rover, as mentioned above, Fermi liquid behaviour is observed only above 
the critical temperature for superconductance anyway. In Section ^, I define 
the non-ladder skeleton functions precisely and prove the above statements 
about E^' uniformly in the temperature in perturbation theory. A nonper- 
turbative proof requires the same extensions of the analysis as mentioned in 
the previous paragraph. 

In pEfl ) an asymmetric model, in which the symmetry E(k) = E(—\&) 
does not hold, was introduced, and a proof was outlined that such models 
are Fermi liquids down to zero temperature. The asymmetry of the Fermi 
surface removes the Cooper instability at zero relative momentum q of the 
Cooper pair, i.e., it implies that the four-point function has no singularity at 
relative momentum q = (which is where the usual Cooper pairing comes 
from). The regularity conditions on the self-energy, which are crucial for 
Fermi liquid behaviour, were not verified in p5|. Doing this is quite a bit 
harder than in the (k — > — k)-symmetric case p^ , ^3| . At zero temperature, 
regularity property (1) has to be replaced by (1') : E € C 2 ~ e because the 
selfenergy is not C 2 at zero temperature (the second derivative grows as a 
power of log/?; for a detailed discussion of these problems see j^]). This 
modified regularity property (1') and (2) were proven in perturbation theory 
for a general class of two-dimensional models with a strictly convex Fermi 
surface, which includes the non-(k — > — k)-symmetric Fermi surfaces, in [21]- 
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p3[ . More precise conditions on the dispersion relation that imply absence 
of the Cooper instability also at nonzero relative momentum q were also 
formulated in [ |22|] . It is not sufficient just to have ak-> — k nonsymmetric 
surface to achieve that; one also needs that the curvature at a point on 
the Fermi surface and at its antipode differ except at finitely many points 
(for details, see Hypothesis (HA 1 ) of 22] and the geometrical discussion in 
Appendix C of ||). 



3 The renormalization group equation 

In this section I derive the continuous RGE for fermionic models. I first 
derive it for the generating function, and then turn to the component form 
which is obtained by expanding the effective action in Wick ordered mono- 
mials of the fields. Let T be a finite set, for a function X \— > F(X) from 
T to any linear space let J r dX F(X) = £r Sxer -^PO' where Er > is 
a constant, let S(X, X') = e^bxx 1 -, an d define the bilinear form (/,<?) = 
Jp dX f(X)g(X). Let A be the finite-dimensional Grassmann algebra gener- 
ated by the generators (ip(X), x(X), rj(X))xer and let S ^ X ) ^ e ^ ne fermionic 
derivative normalized such that g^fy) VK^O = $(X,X'); recall that the 
fermionic derivatives anticommute. 



3.1 The RGE for the generating function 

For t > let Ct be an invertible, antisymmetric linear operator acting on 
functions defined on T, i.e. (C t f)(X) = J r dXC t {X, X')f(X') with 

C t (X',X) = -C t (X,X'). (37) 

Let Ct be continuously differentiable in t; denote ^ t = Ct- Let d[ic t be the 
linear functional (Grassmann Gaussian measure) with characteristic func- 
tion 

J ducM) e M)r = e5 (r? ' Ct v)r . (38) 

The integrals of arbitrary monomials are obtained from this formula by tak- 
ing derivatives with respect to r/. The measure is normalized: J d/ic t (ip) = 1. 

Let V(ip) £ A have no constant part, A € C, and G(0,ip) = XV (ip). The 
effective action at t > is 

0M)=log / <l/ic t (x)e 5(W) . (39) 
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Because the measure is normalized, G(t,ifj) is a well-defined formal power 
series in A. By the nilpotency of the Grassmann variables, 
nomial in A (the degree of which grows with A). Thus Q(t,ip) is analytic in 
A for |A| < A (r). 



Proposition 1 Let 



A 



1 



c, 



2 \5ip' StpJr 



dX / dX' 



C t (X,X') 



5il){X') 



(40) 



Then 



and 



0_ 

at 



£(t,ip) = „A Ct g(o,i>) 



(41) 



(42) 



If Q(0, ip) is an element of the even subalgebra, then for all t > 0, G(t,ip) is 
an element of the even subalgebra, and it satisfies the renormalization group 
equation 



0_ 

dt 



2 V # 



5ip 



(43) 



Proof: For any F(tp) G A, define F(j-) by replacing every factor ip(X) by 
in the polynomial expression for F. Then F(ip) = [F(j-) e^ r, '^' r ] r)= o 



(the derivatives j- also generate a finite-dimensional Grassmann algebra, so 
the expansion for F terminates at some power). Since Grassmann integra- 
tion is a continuous operation, and by Eq. (|38|), 



a 0(o,&) 



ri=0 



For any formal power series f(z) = ^ fk% k , 



■n=o 



/(Ac t ) e 



(44) 



(45) 



so e^> Ct ")r+(^)r = e A Cl e M>)r. Since a c< is bilinear in the derivatives, 
it commutes with all factors that depend only on rj and can be taken out 
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in front in Eq. (p[4|). This implies Eq. ((42|). Since Ac t also commutes with 
A Ct , Eq. @ follows. 

If £/(0, ■0) is an element of the even subalgebra, the same holds for Q(t, ip) 
by Eq. (fl^), since every application of Ac t removes two fields. Thus per- 
forming the derivatives with respect to tp gives Eq. (^). ■ 

3.2 The component RGE in position space 

Let V{ijj) be an element of the even subalgebra. The effective action has the 
expansion 

oo 

g(t,i>) = J2\ r GrM) (46) 

r=l 

with Q r (t, tp) polynomials in the Grassmann algebra. As explained in Section 
|3.1| , this expansion converges for T finite, but the radius of convergence Ao 
depends on T and Ct- For the models discussed in Section 0, this means 
that Ao goes to zero in the limit L — > oo and (3 — > oo. 

Assume that C = lim^oo Ct exists and let D t = C — Ct so that Ct = 
—Df The application will be that C is the covariance of the model and Ct 
is part of it, so that in Q(t,ip), part of the fields have been integrated over. 
Dt is then the covariance of the unintegrated fields. 

I expand the polynomial Q r (t,ip) € A in the basis for the Grassmann 
algebra given by the Wick ordered monomials Q,D t (tj;(Xi) . . . tp(X p )), 

m(r) / m \ 

g r (t, 4>) = J2 I G ™(t 1 20 n Dt n i/>(x k ) , (47) 

m=0 p m \k=l / 

where G mr (t \ X\ : . . . , X m ) is the connected, amputated m-point Green func- 
tion and X_ = (X\ , . . . , X m ) . Details about Wick ordering are provided in 
Appendix ^. A short formula is 

n Dt (rP(X 1 )...^(X p )) = e- AD * 1){X x )...rl>(X p ). (48) 

I use the symbol Q.D t {tp{X\) . . . ip{X p )) rather than : ip{X{) . . . ip{X p ) : to 
indicate clearly with respect to which covariance Wick ordering is done, 
because this will be important. 

The G mr (t | Xi, . . . , X m ) are assumed to be totally antisymmetric, that 
is, for all ir G S mi G mr (t \ X^i) , . . . , ^(m) ) = e(i:)G mr (t | Xi,...,X m ), 
because any part of G that is not antisymmetric would cancel in Eq. (p7|). 
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Application of to Eq. ( |47| ) gives a sum of two terms since two factors 
depend on t. By Eq. (|~ 



= A Ct n Dt (^(X 1 )...ij(X p )). (49) 

When multiplied by G mr (t \ Xi,..., X m ) and integrated over Xi, . . . , X m , 
this gives Ac t G(t,ip). Thus the term linear in Q drops out of Eq. p^ ) by 
Wick ordering with respect to -Dt, and Eq. (|43| ) now reads 



(50) 



m(r) / m \ „ - 

m=0 p m \fc=l / 



where Q r (t,ip) is defined by 



(*$*U^)-gmM> 



(51) 



Being an element of the Grassmann algebra, Q r (t,ip) has the representation 

m(r) / m \ 

Q r (t, ^) = ^ / dX Q mr (t 1 20 ( II ^ Xk ^ ) • ( 52 ) 

To obtain the Q m r(t | X_), one has to rewrite the product of the two Wick 
monomials in Eq. (|5l|). This is done in Appendix |C|. The result is 

Proposition 2 Q m \{t \ X) = 0, and for r > 2, 

Q mr (t\X) = J dK mr j dV J dW (-^detvf\v,W)j 

pi pi 

G rniri (t | Xi,V) G rri2T2 (t \W_,X 2 ), (53) 
where j dK mr stands for the sum 

J dn mr (r 1 ,m 1 ,r 2 ,m 2 ,i) F(r 1 ,m 1 ,r 2 ,m 2 ,i) 

= ^ ^ K mim2 i F(r 1 ,m 1 ,r 2 ,m 2 ,i) (54) 



n>>-2>l (mi,m2,i)6M r ., 

r l~ r " r 2~ r 
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with positive weights K mirn2 i = ("^) (™ 2 ), so that dn mr is a positive measure. 
•M-nram ^ the set of (mi, T7t2, i) suc/t i/iaf i > 1, 1 < mi < m(ri), 1 < mi < 
m(r2), mi + m,2 = m + 2i, anc? mi and mi are even. Xj_ = (Xi, . . . , X mi -i), 
X 2l = (X mi _ i+ i,...,X m ) ) V= {Vx,...,Vi),W= (Wi,...,Wi), anrff = 
(Wj,...,Wi), andvf\]^,W_) is the ixi matrix {vf (V,W)) M = D t (V k ,Wi). 

Comparison of the coefficients gives the component form of the RGE 
d 1 

-Q^G mr (t | Xi, . . . , X m ) = —A m Q mr (t | Xi, . . . , X m ) (55) 
where A m is the antisymmetrization operator 

(A m f)(X 1 ,...,X m ) = ^ Y, e(7r)f(X w(1) ,...,X n{m) ). (56) 

7TG5 

The important feature of Eq. (|53|) is that the determinant of the propagators 
appears in this equation. The Gram bound for this determinant improves 
the combinatorics by a factorial. I now discuss the graphical interpretation 
of the equation and the determinant, to motivate why this improvement can 
be regarded as a 'planarization' of the graphs. 



3.3 The graphical interpretation 

The component form of the RGE has a straighforward graphical interpreta- 
tion. If one associates the vertex drawn in Figure || to G mr (t \X\,..., X m ), 
adopts the convention that the variables occurring on the internal lines of a 
graph are integrated, and writes out the determinant as 

i 

detvf\v,W) = e(n)l[D t (V k ,W v{k) ) (57) 

TTSSi k = l 

the right hand side of Eq. (^) appears as the signed sum over graphs with 
two vertices G min (t) and G m2T2 (t), obtained by joining leg number mi — i+k 
of vertex 1 with leg number i—7r(k) + l of vertex 2, for all k S {1, . . . , i}. The 
graph for ir = id is drawn in Figure ^. The expansion in terms of Feynman 
graphs is generated by iteration of the equivalent integral equation 

t 

G mr (t | X) = G mr (0 | X) + ^A m J dt Q mr (t | X) (58) 

o 
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Figure 1: The vertex corresponding to G mr (t) 



with the initial condition G mr (0 \ X_) = — 5 r iV m (X_), where V m (X) is the 
coefficient of £lc(yj){X\) ■ ■ ■ ip(Xm)) in the original interaction. It is evident 
from Figure || that only connected graphs contribute to this sum. 

Note that the only planar graphs appearing in this sum are from ir(j) = 
j + k mod i for k € {0, . . . , i — 1}, and that for k ^ 0, these permutations 
produce planar graphs only if i = m\ or i = 771,2. For all other permutations, 
the graphs arising are nonplanar. For bosons, the determinant is replaced 
by a permanent, and one can permute the integration variables so that the 
derivative of the permanent (and hence the sum over permutations) gets 
replaced by 

i i\ J] A04, w k ) QjP t {y u Wi). (59) 

k=2 

so that the planar graph drawn in Figure |2| is the only one contributing to the 
right hand side. Thus this factor il distinguishes between the combinatorics 
of the exact theory and a 'planarized' theory, in which i\ is replaced by i 2 
(the second factor i comes from doing the derivative in the determinant, see 
Section |j). The 'planarized' theory does contain more than the sum over 
all planar graphs because of the binomial factors (the antisymmetrization 
operation A m does not change the combinatorics because it contains an 
explicit factor l/m\). In the next section, I bound the determinant by 
const 1 and thereby reduce the combinatorics of the fermionic theory to 
that of the planarized theory. 
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Figure 2: A graph contributing to the right hand side of the RGE 



4 Fermionic sign cancellations 
4.1 The determinant bound 

Eq. ( |53| ) already suggests that a determinant bound similar to the one used 



in Lemma 1 of [20] can be applied to the RGE. Before applying this bound, 
the derivative with respect to t has to be performed, and some factors need to 
be arranged to avoid the factor |A| that appeared in |2(|. The reason it does 
not appear here is that only connected graphs contribute to the effective 
action (whereas the partition function itself was bounded in Lemma 1 of 
[p0| ). In the RGE, this is very easily seen without a reference to graphs. 

Since the determinant is multilinear in the columns of the matrix, the 
derivative with respect to t produces a sum of terms where every column 
gets differentiated. Expanding along each differentiated column gives 

§- t detvP(V,W) =j^j^{-ir i 'b t {V h W v )&etvt X \v®, WP) (60) 

with Vp 1 ' = (Vi, . . . , Vi-x, Vi+i, . . . , Vi) and a similar expression for W( l '\ 
The sign is cancelled by rearranging 

G miri (t \ 2Li,¥_) = {-iy- l G miri (t | X^ l \Vi) 
G m2 r 2 (t\W,X 2 ) = (-iy- l 'G m2r2 (t\ W v ,w} l '\x 2 ). (61) 

Upon renaming of the integration variables, the summand becomes inde- 
pendent of I and I', so the sum gives a factor i 2 . Thus 

Qmr(t | X) = J dK mr i 2 J dV dW D t (V, W) j dYj dZ_ (62) 
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detV? 1] (Y,Z)G miri (t | X 1 ,Y,V)G m2r2 (t \ W,Z,X 2 ). 
Let I • I be the norm [O] 



"mi = max 

pS{l,...,m} 



/m 
l\dX q \F m (X 1 ,...,X„ 
9 = 1 
<3#P 



(63) 



Lemma 1 Assume that 



sup 



det£> t (i_1) (Y,Z) <^i_i(t). 



(64) 



|Qmr(*)|| < / dK mr i 2 A t ^(t) I A I |G mm (i)|| |||G m2J . 2 (t)|||. (65) 



Then 



Proof: Let p G {1, . . . , m}. Without loss of generality, let p < mi — i, so 
that A" p is a component of X 1 (the other case is similar by the symmetry of 
the sum for Q mr {t) in m\ and m 2 ). By Eq. (Bjl), 



HQmr(t)| 



&™ r i 2 Vl(t) sup / [J dX, / dVc/>(V,Xi) X(V) (66) 
x p J q=1 J 



with 



and 



HV,X 1 ) = JdY \G min (t \X lt Y,V)\ 



(67) 
(68) 



X (V) = J dW dZdX 2 D t (V,W) G m2r2 (t | W,Z,X 2 ) 
The bound 

G m2T2 {t | W,1,X 2 ) (69) 



X(V) < sup / 
v J 

gives the result. 



D t (V,W) sup / dZdXj 



From now on I assume that A is a discrete torus and that A* is its 
dual. Let AT be a set with \J\f\ = n elements, let T = A x J\f x {1, 2} and 
r* = A* x M x {1, 2}. J\f can be thought of as the index set containing spin 
and colour indices. The last factor {1, 2} distinguishes between the usual tp 
and as discussed in Section 0. 
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Lemma 2 Let X = (x,o~,j) G T, and for every k E A*, let M(k) be a 
symmetric matrix in A4(n,M.) with eigenvalues m p (k) satisfying \m p (k)\ < 1. 
Let ft : A* — > C, and let 



D t (X,X>) = S fi3 , J dk e^"*') (-iy f t ((-iyk)M a y {(-iyk) . (70) 

A* 

Then D t (X',X) = -D t (X,X'), and 

dk\ft(k)\\ . (71) 



\A» 



Proof: Since A* = —A* and (— 1) J = — (— 1) J ', a change of variables 
A; — > (— l) J k implies 

D t (X,X') = d j+f)3 (-iy [ dk f t (k) M^,(k) e^((-i) f s+(-ip-V). (72) 



The antisymmetry of Dt now follows from the symmetry of M. Let H p {k) 
be the spectral projection to the eigenspace of m p (k), so that M(k) = 
^2 p m p (k)U p (k), and denote the scalar product on the spin space by [•,•] 
so that M{k) aa i = [e a ,M(k) e a >], where the e a are orthonormal. Then 



2tt 



D t (X,X') = J dk ^[a t (X)(k,<p,p), Il p (k)b t (X')(k,<p,p)} 

A*i p 
= (a t (X) , b t (X')) (73) 

with A*i = {k G A* : f t (k) / 0}, and 

at{X){k,<p,p) = (-I)-?' e -^(-i)^ e -^ (i - 3 ) |/t(fc) |i eCT 
bt(X')(k,<p,p) = e ik ^' x ' e™' ft{k)m p (k) \ft(k)\'h a , (74) 

Gram's bound [06|, 



i-1 i 
det0 t (<_1) Cr,#) < H ({at(Y k ),at(Y k )) {b t {Z k ),b t {Z k )))\ (75) 
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and 



\m p (k)\ 2 (e a ,U p (k)e a ) < ^(e CT , IL p (k)e a ) < 1 



(76) 



imply Eq. (0). ■ 

Corollary 1 If ft(k) = L>t(k) and M arJ i{k) = 5 aa i (as in the above many- 
fermion systems) then 



\detvf~ l) {Y_,Z)\ < 



dk 



A* 



A (*0 



i-1 



(77) 



4.2 Power counting for point singularities 



In this section, I show some basic power counting bounds for the Green func- 
tions obtained from the truncation that all marginal or relevant couplings 
are left out. 



Theorem 1 Assume that Dt is of the form Eq. (|70|), that 



dk 



A (AO 



< Aie-*, 



(78) 



and that 



IIAII < A 2 e*. 



(79) 



Let f dk mr denote the measure obtained by restricting the sum in J dn mr to 
mi > 4 and ni2 > 4 and replacing G^ rk by v8 rk \ whenever it appears in 
the sum on the right hand side of Eq. (|5^). Then the solution G mr (t) of Eq. 
(|55j) with initial condition G mr (0) = G^l satisfies 



\\G mr (t)\\ < { 



with "fmr defined recursively as 



7rm-e* ( ^~ 2) ifm>6 
74r(l + *) ifm = 4 
Kl2r ifm = 2 



(80) 



+ 3A2 / dK mr i A-, TmiriTm.2?' > - 

m ' ~ 



(81) 
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Proof: Induction in r, with Eq. ( J80| ) and Eq. (81) as the inductive hy- 
pothesis. The case r = 1 is trivial. Let r > 2, and the statement hold for 
all r' < r. By Lemma |, Eq. (|4|) holds with Ai(t) = A\ _1 e _ *^ _1 ) . Thus 

\lQmr(t)\\ < A 2 f dkmr^Al-h-^WG^m \\G m2r2 (t)|. (82) 



By the inductive hypothesis and mi + ttt-2 — 2i = m, ||| • ||| of the right hand 
side of Eq. fl58|) is bounded by 



||G mr (0)||| + J dkmri 2 A\- 1 A 2 -fm 1 r 1 J m2 r 2 \ f ds e s( f ~ 2) (83) 





To complete the induction step, this has to be bounded by the right hand 
side of Eq. @. If m > 6, f - 2 > f > 1, so 

1 /cfa e s (f - 2 ) < -J— e *(f - 2 ) < A e *(f - 2 ). (84) 

2 y m — 4 m 
o 

If m = 4, | /* ds < £(1 + i). If m = 2, § - 2 = -1, so ± /* ds e s( f ~ 2 ) < 
| = ^. This implies Eq. with 7 mr given by Eq. (|8l]). ■ 

Remark 1 In graphical language, the truncation removes all two-legged and 
all nontrivial fourlegged insertions ('nontrivial' means that the four-legged 
vertices are still there). 



Remark 2 The solution to Eq. ( pi] ) is bounded by the solution to the un- 
truncated recursion 7 m i = v5 m i, and for r > 2, 



7mr = 7 



JS + A— / dK mr i 2 B l 1 7m 1 r 1 7m 2 r 2 (85) 



wii/i ^4 = 3A2 and B = A\. The constants A and B can be scaled out 
of the recursion. If the initial interaction is a four-fermion interaction, 
9m,r = rnjmrA^ 1 B 1 ^^ satisfies the recursion g m ,\ = wb m <±, with w = 4AB, 
and for r > 2, 

r-l 2(s+l) _ / m + 2 fc + 1 _ »\ 
ff m ,r = ^^ ^ ( fc j ( fe Ws ff m+2 fc+2-At,r-s • (86) 

0—1 i.^n u=2 \ / \ 1 / 



S=l fc>0 M=2 
/i even 
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/ do not provide bounds for the solution g mr here; if g mr < const r , then 
the above bounds imply that ^2 r X r G mr (t) is analytic in A. This behaviour 
is suggested by the absence of the factor i\ that would appear for bosons in 
the recursion; see [16 J. 

The propagators Ct and Dt for the renormalization group equation are 
defined using a partition of unity xi +X2 = 1, Xi £ C°°0^o~> [0, 1]) with 

»w-{J (87) 

x'Ax) < for all x E (|, 1), and Hx'illoo < 2 - 

Proposition 3 Let A be d- dimensional, and B be the infinite-volume limit 
of A*. For instance, for A = eL d jUL d , with ^ G 2N, B = R d /^Z d . Let D t 
be given by Eq. flf(j) with \\M(k)\\ < 1, / t (0) = 0, and /or fc ^ 

ft(k) = f(k)xi(e 2t f(k) 2 ), (88) 

uni/i / £ C d+1 (B \ {0},C), satisfying for all \a\ < d + 1 and all \k\ < 1, 
\D a f{k)\ < F d \k\ l - d -\ a \ . Then Eq. © and Eg. © aoZd. 

Proof: Eq. ( ]78|) holds because X\*\{o} \k\ 1 ^ d Xi{k 2 e 2t )dk is a Riemann 
sum approximation to the convergent integral J B \k\ 1 ~ d xi(k 2 e 2t ) ■ The 
infinite-volume analogue of Eq. ((7^) is usually proven by integration by 
parts, using repeatedly 

(x u - x' u )D t (x - x>) = -J ^ d e^-^ k ±-D t (k), (89) 

B 

which implies that \D t (X,X')\ falls off at least as (1 + e~*|x — x'l) - ^ 1 for 
large \x — x'\. On the torus at finite L, one iterates instead the summation 
by parts formula 

D(x - x')(l - e 1 "^')) = fdk e ikix - x,) (t)(k) - D{k - a)) (90) 



which holds for all a £ A*, decomposes A into 2 d parts and chooses a appro- 
priately to get an analogue of Eq. (|89|). This works uniformly in L because 
M(0) = 0. A similar argument is given in more detail in the proof of Lemma 

! ■ 
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Remark 3 Let d > 2, A = eZ d /LZ d and 

2 d 

L £ (k) = 5^(1 — cos(ek v )). (91) 



77ie choices f(k) = L £ (k) 2 , M = 1 (the toy model of ft2C\j), and 

= fe Qsin(M) 2 + (^(^)) 2 ] (92) 



and 

M(k) = f(k) (^i^sin(Ke)+eL £ (k)\ (93) 

(Wilson fermions) satisfy the hypotheses of Proposition^. For all t > 0, the 
infinite-volume limit L — > 00 and the continuum limit e — > 0/ £/te Green 
functions G mr (t) exist and satisfy Eq. (|80[). /n the first case, f(k) — > 
m i/ie second case, M{k)f{k) — > ^ as e — > 0. T/ie Euclidean Dirac matrices 
satisfy the Clifford algebra 7^7^ + 7^7^ = ^z,, and can be chosen hermitian, 
7* = 7^ . is possible to adapt the matrix structure of Lemma [| to satisfy 
the antisymmetry condition on Dt also for this case. 



Remark 4 In ultraviolet renormalizable theories, the signs in the exponents 
of Eq. ( |78| ) and Eq. ( |79] ) are reversed. A power counting theorem similar to 
the infrared power counting Theorem can be proven provided that in the 
initial interaction, at most quartic interactions appear. The statement is 
then that for m>6, the Green function is bounded by const e - *^ -2 -*. 



5 The thermodynamic limit of the many-fermion 
system 

In this section, I apply the RGE to the many-fermion systems defined in 
Section ||| For d = 1, the singularity is pointlike, so the power counting 
bound Theorem [l] applies. For d > 2, the analogue of Theorem || gives 
only weaker bounds, which, e.g., in d = 2 would mean that the four-point 
function is still relevant and the six-point function is marginal. This is not 
the actual behaviour; showing better bounds in d > 2 requires a refinement 
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using the sector technique of [2C] and is deferred to another paper. I also 



show a simple bound for the full Green functions that takes into account the 
sign cancellations. If the coefficients 7 mr given by Eq. ( p5| ) are exponentially 
bounded, this bound implies that the unrenormalized expansion converges in 
a region \X\j3 d < const . I also give a simple proof that the thermodynamic 
limit exists in perturbation theory. 

One can also show bounds on the expansion coefficients for finite n T and 
L. Since this is mainly a tedious repetition of the infinite- volume proofs 
with integrals replaced by Riemann sums (it also requires that /U is chosen 
such that the Fermi surface contains no points of the finite-L momentum 
space lattice), I will content myself with indicating where this is necessary. 

For convenience, I call the limit n T — > oo and L — > oo the thermodynamic 
limit, although the first limit would more aptly be called the time-continuum 
limit. The limit n T — > oo has to be taken first, because I want to apply Eq. 
(||) for operators on a finite-dimensional space only. However, it will turn 
out that for T > 0, the order of the two limits does not matter. Since I want 
to take the limit n T — > oo, I can assume that n T > 2/3(eo + Pmax)- Thus 
Lemma ^ applies, with Eq = et, for all t > 0. 

5.1 The component RGE in Fourier space 

Let A = T x L and A* = M„ T x L*, as given in Sections [O] and ^ Let 
r = A x {-1,1} x {1,2} and T* = A* x {-1,1} x {1,2}. The Fourier 
transforms of the G mr (t) are, with (P = (Pi, . . . , P m ), 

/m 
\{ dx k e-^ + - + ^ x ^G mr (t | X). (94) 
J i 



k=l 



For K = (k,a,j) G T* let ~ K = (-k,a,3 - j), and let 5 r *(K,K') = 
8A*(k, k')5 aa t5jji . Assume that the Fourier transform of the propagator D t 
is of the form 

Dt(K, K') = 8 V * (K, ~ K') D t (K') (95) 

where 

D t (K) = (-iyD t ((-iyk). (96) 

This combination of signs implies that Dt(X,Y) = —Dt(Y,X). The prop- 
agator Dt for the many-fermion system is given in Eq. (|103|). The Fourier 



29 



transform of Q mr (t) is 

Q mr (t | P) = j dn mr Hi ( dK x ... dKi (-Dt(Ki)) f[ D t {K 3 ) 

J J 3=2 

G miri (t | P {1 \K) G m2r2 (t \~K,pW) (97) 

withpW = (Px,...,P mi ^),P^ = (P mi _ i+1 ,...,P mi+m2 ),K = (K 1 ,...,K i ), 
and ~ K_ = (~ . . . , ~ K\). Here the arguments of G m2r2 (t) have been 
permuted and relabelled such that the determinant is transformed into a 
sum with the same sign for all terms (see Appendix Q); this gives the extra 
factor i\. 

By translation invariance in space x and time r, 

G mr (t | Pi,..., P m ) =5 A .(pi + ...+ Pm , 0)I mr {t I Pi,...,P m ) (98) 

with a totally antisymmetric function I mr (t | Pi, ... , P m ) of (Pi, . . . , P m ) € 
T* m that satisfies 

m 

J2^ P Jmr(t\Pl,...,P m )=0. (99) 
M=l 

A priori, Eq. ( |98| ) implies only the existence of a function I m r(t | Pi, ... , P m ), 
defined only for those Pi, ... , P m for which pi + . . . +p m = 0. However, since 
H = {pi, . . . ,p m : p\ + . . . +p m = 0} is a linear subspace of (A*) m , one can 
simply extend / to a function on all space by defining I = IoTl^ with TIh the 
projection to the subspace H. Since Hff is symmetric in all its arguments, / 
is totally antisymmetric, and Eq. (|^) holds because (1, . . . , 1) _L H (in Eq. 
(|99|), V Pm is a difference operator which becomes the gradient in the limit 
where the momenta become continuous). 

The product of the two <5a* in Eq. ( |97|) can be combined to cancel the 
5 a* in the relation between G and /, and to remove the integration over k\. 
Thus the RGE in Fourier space is J^Pmr^ I P) = \^mQmr{t I P), with 

Qmr(t | P) = J dK mr Hi J dK 2 . . . dK % £ (~A(^l)) 

i 

H D t (K s ) I min (t | P (1) ,iQ Im 2 r 2 {t |~ K,P( 2 ))(100) 

where i^i = (fcj, 0"i,ii) and /ci is fixed as A;i = — (&2 + . . . + fcj) + p\ + . . . + 
p mi _i+i. I use the same symbol for the Q in position and in momentum 
space since it will be always clear from the context which one is meant. 
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Remark 5 The equivalent integral equation is 

t 

I mr (t | P) = I mr (0 | P) + ^A m J ds Q mr (s | P). (101) 

o 

The sum J d,K mr contains the sum over r\ > 1 and r2 > 1, with the restric- 
tion n + r2 = r. XTms only ri < r and r2 < r occur in this sum. Therefore 
Eq. (|101| ), together with an initial condition 7 mr (0 | P), uniquely determines 
the family of functions I mr (t | P). Iteration of Eq. ( |101| ) generates the usual 
perturbation expansion. 



5.2 Bounds on the finite-volume propagator 

The propagator c for the many-fermion system is given in Eq. (|l7|). If k 
is such that E(k) = 0, then c becomes of order /3 for small In the 
temperature zero limit, j3 — > oo, this becomes a singularity. This is the 
reason why renormalization is necessary. The renormalization group flow is 
parametrized by t > 0, where 

e t = eoe"* (102) 
is a decreasing energy scale. The fixed energy scale eo was specified in 



Section 2.3. The limit of interest is i — ^ oo. 

The uncutoff propagator for the many-fermion system is given in Eq. 
Thus, for k = (w,k) € A* let 

A(fc) = ^ _^ (k) Xi (V 2 |i£ - ^(k)| 2 ) (103) 

where Q is defined in Eq. (0), Xi is given in Eq. (|87|), and define Ct(fc) 
similarly, with xi replaced by \2 = 1 — Xl- Then Ct(k) + D t (k) = (iu; — 
-E(k)) -1 is independent of t. The functions Ct and A define operators 
I) t {K,K') and C t (K,K r ) on the functions on T* by Eq. (p§) and Eq. 



Denote A = -§jDt, and let 11 (A) = 1 if the event A is true and otherwise. 
Proposition 4 Df is a C°° function of t that vanishes identically if t > 



supp A C {k € A* : |iu) - £?(k)| < e t } 

suppA C {fcG A* : -e t < |icD — £7(k)| <e t }. (104) 
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Moreover 



D t (k) <4e t _1 l(\iu-E(k)\ <e t ) < 2/3 l(\iQ - E(k)\ < e t ) , (105) 



D t (k) < § l(|i£-£(k)| <e t ) , 



dk 



D t (k) 



< 4Vi and / dk 

A* A' 

where Vi is the constant in Eq. ([l9|). 



AW 



< Vi log ^ 



2 ' 



(106) 



Proof: Xi( x ) = if x > 1, so Dt(k) ^ implies |icD — E(k)\ < e^ < 1. By 
Lemma [|, this implies |w| < fe^. Since |w| > jj, Dt ^ ® onr y f° r * < 1°§ 2 • 
Since |i£D 

t-derivative gives \%Dt{k)\ = 2e~[ 2 \iQ 



E{k)\ > | Re cD| > the stated properties of Dt follow. The 

E(k)\ \ X [^ 2 \iQ - E(k)\ 2 )\. Since 



Xi(x) = for x (±,1), Z? t (fc) / implies \e t < fi£D - £(k)| < e t which 
implies Eq. ( |105| ) and Eq. ( |104| ). Eq. ( [L06| ) follows from these inequalities 

by Lemma ||, by Eq. (|l^) , and by D t = J t 



log(/3e /2) 



(is ZX 



Remark 6 T/ie bounds Eq. (|l06j) are crude because the restriction \E(k)\ < 
e+ was replaced by \E(k)\ < 2 w/ien Eg. @ was applied. To get a better 
bound, one has to require that no points of the finite-volume lattice in mo- 
mentum space is on the Fermi surface S. Because only then J L dk 1 (|-E(k)| < et) < 
const et holds uniformly in L. 



5.3 Existence of the thermodynamic limit in perturbation 
theory 

The proof of existence of the thermodynamic limit will proceed inductively 
in r, because of the recursive structure of the RGE mentioned in Remark 
||. It will be an application of the dominated convergence theorem to Eq. 
(|101|). To this end, it is necessary to make the integration region independent 



of n T and L. Although I mr (t), given by Eq. ( |101|) , appears evaluated at 
P = (Pi, . . . , P m ) E T m on the RHS of Eq. ( [TOl]) only, the integral defines 
the t-derivative of a function defined on T^" 2 , where = M(J3) x B x 
{1, —1} x {1,2} with B = G* the first Brillouin zone of the infinite lattice, 
and 

M(fi) = {uj n = ^(2n + 1) : n € Z} (107) 
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the set of Matsubara frequencies in the limit n T — > oo. For a bounded 
function F m : r^ m -» C, let 

I -^m I o — SUp | F m (P)| (108) 

p<=r* m 

no 



Up to now, the dependence of I mr (t) on (n T ,L) was not denoted explic- 

lim = 

(n T ,Z/)— ► oo 



itly. I now put it in a superscript and write Lmr' L \ Let lim 



lim^oo lim nr 

— »oo • 

Lemma 3 Let (Imr' L \o)) m ,r,n T ,L be a family of bounded functions such that 
Imr' L \o) = if m > 2r + 2, ]xmr nr Imr' L \o) = ^mr(O) exists and is a 

bounded function on r^ m , and Lmr' L \o) < Kml- Let (Imr' L \t))m,r,n T ,L be 

the solution to Eq. ( |101| ). Then, for allm andr, Imr' L \t) = ifm > 2r + 2, 
and 

^Lt T ' L) (t) = for all t > log |, (109) 
there are bounded functions L mr (t) : T^" 1 — » C smc/i i/iaf 

lim 4"- L )(t)=/ mr (t). (110) 

(n T ,L)— ioo 

Let P mr 6e i/ie polynomials defined recursively as 

P mr (x) = 6 1-r |I mr (0)| + eg 1 J dn mr i i\ x l ~ 1 P miri (x)P m2r2 (x) (111) 

(in particular, the coefficients of P mr are independent of (3), then for all 
n T ,L,(3,t 

lt' L) (t) „ < (eoPy-'Pmr (m log Sf) . (112) 







Proof: Induction in r, with the statement of the Lemma as the inductive 
hypothesis. Let r = 1. Since r\ > 1 and r2 > 1, the right hand side of the 
equation is zero, so 1^ (t, P) = ^'^(0,?) for all i. Thus the statement 
follows from the hypotheses on Lmr' L \o). 

Let r > 2, and the statement hold for all r' < r. Let m > 2r + 2. Then 
J^' L) (0) = 0. The five-tuple (mi , n, iri2, r 2 , i) contributes to the right hand 
side only if n + r2 = r, mi + m 2 = m + 2i, i > 1, and by the inductive 
hypothesis, only if mi < 2n + 2 and wi2 < 2r 2 + 2. Thus L^> L \t) can be 
nonzero only if m = m\ + m 2 — 2i < mi + m2 — 2 < 2ri + 2 + 2r 2 + 2 — 2 = 
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2r + 2. The integral appearing on the right hand side of Eq. ( |101| ) is a 
Riemann sum approximation to an integral over the (t, n T , L)-independent 
region s E [0, oo), kj £ M(/3) x £?. By the inductive hypothesis, the factors 
Imk^ have a limit Eq. ( p.!0|) satisfying Eq. ( |112j ). Since for all i, Z) t is 
a bounded C 2 function, the same holds for the propagators (boundedness 
holds because P < oo). Thus the integrand converges pointwise, and it 
suffices to show that it is bounded by an integrable function to get Eq. ( |110| ) 
and Eq. ( |112|) by an application of the dominated convergence theorem. Let 



a = 4Vi log and let g be the function on [0, oo) x (M(/3) x B) 1 1 given 



2 

by 



r 2 -l 



g{s,k) = — s 1 (s < log ^a) y dn mr i i\P miri {a)p r ^ l P m2ri {a)p 
1 f 

Yl / dsj 1 (|cjj| < f e e - ^') 11 (|S(kj-| < 2e e - ^) . (113) 

The integrand on the right hand side of Eq. (|101|) is bounded by g by Proposi- 
tion Lemma 0, and the inductive hypothesis Eq. (|112j ). Because g vanishes 
identically for s > log^, it is integrable. By Eq. (p.05|) and Eq. (|106|). 



t 

ds I g(k)dk 2 ...dh< p r ~ l I dn mr i i\ a l ~ l P miri (a)P m2r2 {a). (114) 



o 



Thus Eq. ( pq ) holds, and Eq. ([112]) holds with P mr given by Eq. (|ITlD (the 
factor 6 1_r comes from the assumption that Pe$ > 6). Since the right hand 
side of Eq. (|ToI| ) vanishes for t > log ^& Eq. (|i%9|) holds. ■ 



Remark 7 /m r r ' (0) is not the initial interaction because at t = 0, the 
propagator Cq ^ 0. The Imr' L \o) are obtained from the original interaction 
by Wick ordering with respect to C and integrating over all fields with co- 
variance Co- The existence of the limit (n T , L) — > oo of the Imr' L \o) is not 
obvious because in the limit n T — > oo, the absolute value of the propagator 
is not summable. In perturbation theory, this is no serious problem, and 
the Imr' L \o) are controlled by a similar inductive proof as the above. An 
alternative proof is in Appendix D of 
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5.4 The many-fermion system in the thermodynamic limit 

In this section, let the many-fermion system satisfy the hypotheses stated 
in Section ^3] with ko > d. 

In the thermodynamic limit, the spatial part p of momentum becomes 
continuous, p G B, and the set of Matsubara frequencies becomes M(/?), 



given in Eq. (107). It is convenient to take (po)P) £ R x B and put all the 
/3-dependence into the integrand. To this end, I define the step function 
top : R -> M(/3) by 

^(p ) = -(2n + l) ifp Q G(— n,— (n + 1)]. (115) 

For any continuous and integrable function /, 

J ^ /( W/9 (po)) = ^ E /(«)■ ( 116 ) 



Moreover, 



IT 7T 

sup \top(p ) -po\ = ~s and inf \iop(p )\ = -5, (117) 
p eR P P06M p 



\up(po) I > and top(-po) = — iop{po) holds Lebesgue-almost everywhere, 
so that in integrals like Eq. ( |116| ), top can be treated as an antisymmetric 
function. With this, the propagator now reads Dt{p) = Ct(iop(po), E(p)) 
with 

C t (x,y) = -J— Xl (el 2 {x 2 + y 2 )) (118) 
ix — y 

(and e t = e e~*). 

In infinite volume, the restriction on the spatial part p of momentum 
improves the bounds on the propagator over the finite-volume ones given 
above. 

Lemma 4 D t is bounded, C°° in t and C ko in p. If t > log^ 2 -, then 
Dt(p) = for allp GRxB. For all multiindices a € Nq with \a\ < ko, there 
is a constant B a > such that 

D a £) t {p)\ < B^ 1 '^ l(\iup(p ) - E(p)\ < e t ) 

< B^ 1 ^ t{\up(p )\ <e t ) t(\E(p)\ <e t ) . (119) 
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S = 4. Fort < log £a 



2?r 



(120) 



and 



(27r) d + 1 



l-|a| 



/ 



(27T) d + 1 



D t (p) 



< 8Jie t . 



(121) 



(122) 



Proof: For .Dt(p) to be nonzero, ^(po)! < et must hold. By Eq. ( |117| ), 

D t = Ofort > log^a. Since = -2e t - 2 (^(po)+£(p)) xi(^ 2 (^(Po) 2 + 

£(p) 2 )) and Wx'iW^ < 2, |A(P)| < | l(|^(p )-^(p)|<e t ) . Thus 
-Bo = 4. Derivatives with respect to p can act on -E'(p) or on Xii m which 
case they produce factors bounded by e^~ 1 \VE(p)\, so there is B a such that 
Eq. ( |119| ) holds. Inserting Eq. ( |119| ) into the integral in Eq. ( |120D gives 



dpp 
2vr 



D a D t (p) 



< B a e t 



-l-\a\ 



t{\E(p)\ <€ t ) Mg 



(123) 



with 



{n G Z : |2n+ 1| < — } 

7T 



2 

< -e t , 

IT 



which proves Eq. (|120|). This gives for the integral in Eq. 



d d+1 



P 



(2vr 



D a D t (p) 



< B a e t 



-H 



P 



(2vr) c 



l(|#(p)|<e t ) 



(124) 



(125) 



IRxB B 
For t > 0, et < eo, so by a change of coordinates in the integral, 

j P l(|^(p)| <e t ) < j dp f d9\J(p,0)\ = 2,he t (126) 



which implies Eq. (|12l| ). Eq 



follows by integration over t. 
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5.5 A bound in position space 

In this section, I prove a bound for the unrenormalized Green functions 
that motivates the statement that the unrenormalized expansion converges 
for \\\(3 d < const (if the solution to Eq. (|86|) is exponentially bounded, it 
implies this convergence). 

Lemma 5 For many-fermion systems with E G C k ° , where ko > d, there is 
A2 > such that for all t > and all (5 



All < A 2 e* d . (127) 

1/2 

0/2 



Proof: By definition of ||| • |||, |||-Dt||| < 4 f^, 2 dxo f^dx. |D((xo,x)| with 



(x ,x) = /tf / ^e ixo ^+ ik -%(up(k ),E(^ (128) 



a^(fc ) 2 +£(k) 2 



By Eq. O, C^ko), E(k)) = -£(i^(*o) + S(k))xi( 
Claim: There is a constant AT, depending on A:q and -E, such that for all 

\ti> t (xo,x)\<Net — — 1 : (129) 

(l + e t |x|) fc o (l + et min{|x |,f - |x |}) 2 

The lemma follows from Eq. (|129|) by 



and 



/ ^ (l+ e l|x|)*° = ^ I ^ (1 + lel) " fc0 (130) 





2 OO 

1 i f du , 

dx 5 < / T, v»- (131) 

(l + e t min{|x |,f - |x |}) 2 " 7(l + n) 2 



Eq. (129) is proven by the standard integration by parts method. Since the 
/co-dependence is via the step function up, one has to use summation by 
parts in the form 

B t (x ,x)(l-e 2 -^y = | fa J I ^e i ^( fc o) +ik - x (ALCi)(^(fc ),ii;(k)) 

R B 

(132) 
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where (A a /)(&o) = f{ko) ~~ f(ko — a)- Taking q < 2 and using a Taylor 
expansion of Ct(x,y) to order g in x, one sees that A^Ct is still C fc ° in k 

and that for any multiindex a with |a| < k$, 



(^rAtC t (^(k ),E(k)) 



<Ni - l(|^(fco)-£:(k)|<e t ) 

e t 



By Eq. fll2j) and Eq. (|||), this implies 



x' ! ! ! -eT iX0 ) 2 D t (x) 



For |e| < 7T, | sin£| > £ min{|£|,7r - so 



1-e" 



> 



2vr 
sm — x 

/5 



> -min{|x |, § - |x |}, 



and thus 



6[ a|+ V|(min{|* |,f -ko|}) 9 
This implies Eq. ( |129| ), 



< N 2 e t . 



(133) 
(134) 



(135) 
(136) 



Remark 8 This proof changes in finite volume: one needs the requirement 
that no point on the Fermi surface S is a point of the momentum space 
lattice determined by the sidelength L to get a bound that is uniform in L. 
This can be achieved by choosing the chemical potential fj, appropriately. 



Theorem 2 For the many-fermion model in d > 1, and with the initial 
condition |||G mr .(0)||| = 7m', 



II GW(t) I <7mr(eoP) (r - 1)d 
with the j mr given by Eq. flB5|) ; where A = -^^2 and B = 8Jieo- 



(137) 



Proof: By Eq. (ggp, Corollary |, and Lemma g | detX> t (i 1} | < (SJjet)*- 1 , 
so Ai < 8Jieo- By Lemma [5], \\Dt\\ < A2e* d . The proof is by induction in 
r, with the inductive hypothesis 



l|Gmr(t)|| < 



td(r-l) 



(138) 
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The statement is trivial for r = 1. Let r > 2, and Eq. Q13SD hold for all 
r' < r. By Lemma [j] and Eq. (|106| ), and since r\ + r 2 = r, 

\lQmr(t)\\ < A 2 e td ^ J d Kmr i 2 A*" 1 7m 1 n7m a n«- ( 139 ) 

Since r > 2, r - 1 > §, so Jj dt e^" 1 ) < ^e^" 1 ), and Eq. Q) follows 
by integration and by r > ^. Because = if i > log(/3eo), Eq. (|137| ) 
holds. ■ 



Remark 9 The bounds here are rather crude because et < eo was used. 
However, these are bounds for the full Green functions, not for truncations. 
Improving them requires renormalization and, ford > 2, the sector technique 

of El- 



Proposition 5 The hypotheses of Theorem || are satisfied for the many- 
fermion model in d = 1 in the thermodynamic limit. 

Proof: Eq. (ff^) follows from Eq. (|121| ). Eq. (^) follows from Lemma |. 



Proposition 6 Let d > 2. If the sum J dK mr is truncated to mi > 2(d+ 1), 
and m,2 > 2(d + 1), then the G mr defined by this truncation satisfy for all 
m > 2d + 2 

|G mr (i)l <lmr^- {d+1)) (140) 

If these bounds were sharp, all m-point functions up to m = 2d would be 
relevant in the RG sense, and the (2d + 2)-point function would be marginal. 
This is not really the case; see |S6|/ . It is, however, remarkable that such a 
simple power counting ansatz goes through at all in d > 2. 



6 Regularity of the selfenergy 



In this section, I study the regularity question for the skeleton self-energy, 
which is very nontrivial even in perturbation theory. I show that all notions 
23] arise in a natural way in the RGE and prove that 

This verifies the regularity 



introduced in [21, 22, 
the skeleton selfenergy is twice differentiable. 
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criterion (1) for Fermi liquids in perturbation theory. The proof given here 
is a considerable simplification of the ones in |2l] , . 

The main technique in doing the regularity proofs, and in showing that 
the generalized ladder graphs give the only contribution to the four-point 
function that is not bounded uniformly in /3, is the volume improvement 
technique invented in [ffifj . Here, I show that the overlapping loops, and 
all related concepts developed in |]2l|| -p^1, appear naturally in the Wick- 
ordered continuous RGE. In particular, overlapping loops always appear in 
the skeleton two-point function, and the only nonoverlapping part of the 
skeleton four-point function is mi = m<i = 4, which is precisely the ladder 



part of the four-point function. Moreover, the double overlaps of [23| arise 



in a natural way when the integral equation Eq. (101) is iterated. The 
main reason why all these effects are seen so easily is Wick ordering. The 
RGE without Wick ordering has too little structure to make these effects 
explicit in a convenient way. For instance, the volume improvement effect 
from overlapping loops is a two- loop effect, whereas the non- Wick-ordered 
RGE is a one-loop equation. 

6.1 Volume-improved bounds 

The indicator functions restrict the spatial support of the propagator to 
regions 

TZ(e t ) = {p e B : \E(p)\ < e t }. (141) 

The volume of the intersection of 1Z and its translates occurs in the bounds 
for the integrals in the RGE. Good bounds for these volumes are the key to 
the analysis of these systems. 



Lemma 6 Let ko > 2, ir(p, 6) be as defined in Section \2.$ , and 

W(e)=sup max / dQ x \ d0 2 t (\E(vnr(0, 6i) + «i7r(0, Q 2 ) + q)| < e) 
q eB^e{±i} J J 
S d-i S d-i 

(142) 

There is a constant Qy > 1 such that for all < e < eo 

M,)<^{; +|log£| :^>3. <« 3 » 

Proof: This is Theorem 1.2 of The constant Qv depends on the 

curvature of the Fermi surface, hence on eo- ■ 
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I now prepare for the regularity proofs by using the volume improvement 
of Lemma |6| to bound the right hand side of the RGE. This is possible 
because the graph drawn in Figure ^ is overlapping according to the graph 
classification of [21], if % > 3. 

D a Q mr (t\P) = j dK rnr {m l ,r l ,m2,r 2l i)ii\X a ,i(t\ P) (144) 



where, after a change of variables kj — > (— 1) j k 



3i 



*e{-i,iy ie{i,2y a RxB 

i 

H Dtikj) D a U miri (t | PP-),K) D a *I m2r2 (t \~K,PW) (145) 

3=2 

with Kj = (kj,(Tj,lj), the sum over a running over all triples (00,01,02) 
with ao + ai + 02 = a, and k\ = - 2} =2 (— l)' 1 '^" +Pi + - ■■ +Pm 1 -i+i- %a,i 
depends also on {m\,r\, m 2 , r 2 , m, r). Application of | • | gives 

\XaMo < 4 'E^ra \D ai I min (t)\a \D a H m2r2 (t)\ y ao>i (t) (146) 

a 

with 

/i i 
\\ \D t (kj)\dkj\D a D t (T V jkj + Q)\ (147) 
7 . , 



3=2 3=2 



Lemma 7 y a ,i(t) = for t > log and 

y a ,i(t) < (SJi)* -1 ^ eT 2 " |a| . (148) 

//»>3, 

y a ,i(t) < (W^BaK^l + t) er 1Ha| (149) 
<7?uen in (|158[) . 
Proo/: By Eq. (fill), |D a A(p)| < B«er Ha| , so 



i-1 



!M<) < B a e t 1 |a| y |A(fe)|rffc] , (150) 
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thus Eq. (pgp holds by Eq. (|i"2l|) . For i > 3, 



i-3 



3Wt)<[ / |A(fc)|dfc| sup max FofoQ.vi,^) (151) 

\ J I QeM(/3)xB"i^2G{-l,l} 

\MxB / 

with 

Y a (t,Q,v 1: v 2 ) = JdhJ dk 2 \D t (k 1 )\ \D t {k 2 )\ \D a D^v^ + v 2 k 2 + Q)\. 

(152) 

In Y a (t), I use D t = for t > log to write D t (kj) = - J^ eo ^ dtjD tj (kj). 

By Eq. ( |TT9|) , \D a D t (p)\ < B^ 1 ^ l(|£(p)| < e t ) . Inserting this and 
doing the integrals over (A;i)o and (&2)o> I get by Eq. Q120 ) 

Log(/3e ) log(/3e ) 

r a (t,Q,m,«2) < 8 2 /3 Q e~ Ha| /" dti j dt 2 V(t,h,t 2 ) (153) 



with 



V(t,ti,t2)= J dki y dk 2 Id^Cuiki + uaka + Q)! <e e _t ) 

^(eoe-*!) 7£(e e-*2) 

(154) 

With the coordinates (p, 0) defined in Eq. (pl|), 



e e "i e e 



-'2 



v(t,ti,* 2 ) = j dpl j d P2 J j de 2 j( Pl ,e 1 )j(p 2 ,e 2 ) 

— eoe _t i — eoe — * 2 

l{\E(v 1 n(p u e 1 )+v 2 w(p 2 ,6 2 ) + Q)\<e e- t ) . (155) 

Since < eoe~ tj < eoe - *, Eq. ( pi|) implies ir(pi, 0i)+t> 2 7r(p 2 , # 2 )+Q) 
-S(«iir(O,0i) + W2ir(O,02) + Q)| < ^l^ljeoe-*, with \E\ Y = \VE\ . Using 
J(/?j, Oj) < Jo and doing the p-integrals, I get 

V(t,h,t 2 ) < (2e J ) 2 e-* 1 - te W ((1 + ^f^oe"*) (156) 

with W given by Eq. ( |142| ) . The integrals over t\ and t 2 give 

Y a (t,Q, Vl ,v 2 ) < 4 4 J B Q J 2 eJ-^W ((1 + i g i )e e-*) . (157) 
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This implies Eq. ( |149| ) with 



Kl*> = 2(^) 2 Qv(l + ^) (l + log(l + ») + llogeol) . (158) 

Jl On V 90 J 



This gives the following 
Lemma 8 For all i, 

\X a ,i(t)\ < 4H8Jir 1 ^^^ T | J D ai W 1 (t)| | J D a2 W 2 (t)lo 

a 

B ao ei~ 2 ~ la ° l - (159) 

For i > 3 

a 

B^Wej- 1 - 100 '/^^*) */ d = 2 (160) 
* \l ifd>3. ( J 

Eq. ( |160| ) is the implementation of the volume improvement bound from 
overlapping loops in the continuous RGE setting. In the next section I use 
it to prove regularity properties of the selfenergy. 

6.2 The ladder four-point function and the skeleton selfen- 
ergy 

Definition 2 The RGE for the skeleton functions I mr (t)is given by Eq. 

), but with f dn mr replaced by J dh mr , where in the latter the sums over 
mi and rri2 start at four instead of two. The function Izrif) is the skeleton 
selfenergy of the model. 

This truncation of the sum prevents any two-legged insertions from oc- 
curring in the graphical expansion. Thus only graphs without two-legged 
proper subgraphs contribute to the skeleton functions. Since all effective ver- 
tices occurring in our model have an even number of legs, the skeleton self- 
energy S r = /2,r is indeed given by the sum over all two-particle-irreducible 
graphs (this is proven in Proposition 2.7 of [^31 ) and the skeleton four-point 
function is the one-particle-irreducible part of the four-point function (this 
is proven in Remark 2.23 of pi]]). 
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Figure 3: The graph corresponding to Q 4 



r,2 



To give a precise meaning to the statement that the ladder resummation 
takes into account the most singular contributions to the four-point function, 
I split the skeleton four-point function into two pieces, as follows. In Q 4r , 
m\ + m-2 = 4 + 2i, so i > \{rn\ + m 2 — 4) > 2 since the skeleton condition 
has removed m\ = 2 and 771,2 = 2. Let Q 4jr ,2 be the i = 2 term in this sum; 
it corresponds to the 'bubble' graph drawn in Figure [| 

More explicitly, Qa^,2 is given by 



1 



Q&,rfl\t I Pli - ■ ■ 1 Pi) — 



144 E E I dk D t (K) D t (K') 
h, ri (t I Pi,P 2 ,K',K) J 4>ra (t |~ K,~K',P 3 ,P 4 j\ (161) 
with K = (k,ii,ai) and if' = (p% + P2 — k, 12,(72)- Let 

Q 4 ,r,>3(i I E) = QaA* I P) - QiM* I £) (162) 

be the contribution from all terms where at least i > 3 lines connect the two 
vertices in Q 4r . Correspondingly, let 



/ 4 , r (i) = B r (t) + U r (t) 
where -B 4jr (£) and U^ r (t) are defined by 



0_ 

Of 



B r (t) = ^A 4 Q V)2 (i), ^U r (t) = ^A 4 Q V) > 3 (i) 



(163) 



(164) 



with the initial condition B r (0) + U r (0) = / 4)r (0) (where it is understood 
that one of the two summands on the left hand side is set to zero, see below). 
Note that because I 4ir occurs on the right hand side of both equations in 
Eq. (|164|), Eq. (|164j) is a coupled system of differential equations. 
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Definition 3 The function B^ obtained from the skeleton RGE by the 
truncation 

^U r (t) = 0, U r (0) = (165) 

is the ladder skeleton four-point function. The function Ur obtained by 
the truncation 

^B r (t) = 0, B r (0) = (166) 

is the non-ladder skeleton four-point function. The functions Imr> obtained 
with this truncation, (i.e., where m\ = 4 and ni2 = 4 are left out in the sum 
J dR^ r ) are the non-ladder skeleton Green functions. 

The motivation for the split is that in Q4 ir ,>3, the number of internal lines 
of the graph in Figure ^ is i > 3, so the volume improvement of Lemmas ^ 
and H improves the power counting. 

The constants in the following theorems are independent of (3. 

Theorem 3 There is a constant L2 such that 

Bf\t) KWiUl + t^-^Wllog^r 1 (167) 
— ■ 

The series ^ B r (t) converges uniformly in t and JP if | A log /3eo| < I^ -1 - 



Proof: This follows immediately by induction from Eq. ( 161 ) by use of 
Eq. (pp with i = 2 and a = 0. ■ 



Using one-loop volume bounds one can show that both the particle- 
particle and the particle-hole ladder are bounded uniformly in t and (3 for 
Q 7^ 0, so that the only singularity in the four-point function can arise at zero 
momentum (for a discussion of this, see, e.g., [§4|]). In the particle-particle 
ladder, this singularity is really there, and it implies that Fermi liquid be- 
haviour occurs only above a critical temperature: the log (3 in Definition [l] 
is the logarithm occurring in Theorem |3|. 

The next theorem states that the non-ladder skeleton four-point function 
is bounded and that consequently, the non-ladder skeleton selfenergy is C 1 
uniformly in (3. This shows that indeed, only the ladder four-point function 
produces a nonuniformity in (3 (this motivates the alternative criterion for 
Fermi liquid behaviour in Section ^1]). The second derivative, however, is 
only bounded by a power of log/3 in d = 2; this motivates why at zero 
temperature, the self-energy is only required to be C 2 ~ 5 for some 5 > 0. 
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Theorem 4 For all r > 1, the non-ladder skeleton functions I^Pit) con- 
verge for t — > oo to a C 2 function. There are constants L^^ and L^ r , 
independent of 13, such that 



D«~A N r \t) 



< L 



3,r 



and 



D a i ( H\t) 



1 if \a\ < 1 

(log/3eo) 2 if \a\ = 2 and d = 2 
log/3eo if \a\ = 2 and d>3. 

1 if \a\ = 

< L 4>r { log/3e if \a\ = 1 
/5e if \a\ = 2. 



(168) 



(169) 



If the ladder four-point function is left in the RGE, its logarithmic growth 
(Theorem |3|) shows up in all other Green functions, and in the selfenergy: 

Theorem 5 The skeleton functions Ii T and I^ r converge for t — > oo and 
satisfy 

D a I mtr (t)\ o < L 5imjr (log/3e ) r (170) 
for m = 2, \a\ < 2, and m = 4, \a\ = 0. 
Theorems § and || are proven in the next section. 

Remark 10 To prove bounds with a good (3 behaviour for m > 6 requires the 
use of different norms, where part of the momenta are integrated jl^/. It is 
easy to see that form > 6, the connected m-point functions have singularities 
and thus are not uniformly bounded functions of momentum; for instance 
the second order six-point function is given by C(P), which is 0{(3) if p is 
on the Fermi surface and if to = n/ '(3. 



Remark 11 A bound const r for the constants in Theorem would suf- 
fice to show the first requirement of the Fermi liquid behaviour defined in 
Definition [J, namely the convergence of the perturbation expansion for the 
skeleton functions for | A log /3| < const . The proofs given in this section 
do not imply this bound because in the momentum space equation, the fac- 
torial remains. It may, however, be possible to give such bounds in d = 2 by 
combining the sector technique of jj^j with the determinant bound. 
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6.3 The regularity proofs 

At positive temperature, the frequencies are still discrete, whereas the spatial 
part of momentum is a continuous variable. Whenever derivatives with 
respect to po are written below, they are understood as a difference operation 
2^(/(Po + j?) ~ f(Po))- The RGE actually defines the Green functions 
for (almost all) real values of po, not just the discrete set of Matsubara 
frequencies, so the effect of such a difference can be bounded by Taylor 
expansion. This changes at most constants, so I shall not write this out 
explicitly in the proofs. 

The following theorem implies Theorem |i| about the non-ladder skeleton 
selfenergy and four-point function. Note that all bounds in this Theorem 
are independent of (3. 

Theorem 6 Let a be a multiindex with \a\ < 2. For all m,r, let L mr (0) be 
C 2 in (pi, . . . , Pm), and assume that there are Kml > such that |D a / mr (0)| < 
Kml, with Kml = for m > 2r + 2, and = 5 m 4V, where v > 0. Let 
Lmr {t) be the Green functions generated by the non-ladder skeleton RGE, 
as given in Definition [| ; with initial values L mr (0). Let K m \ = and for 



where Mi = 24QBK^ , with B = max a B a . Then K mr = if m > 2r + 2, 
and for all t > and all m, r 



r > 2, let 




(171) 



if m > 6 
if m = 4 
if m = 2 



t 



D 



mr 



2-|a| n +t \S d , 2 



(172) 



Moreover, for m > 6 



rrivpit) <K mr 4 



m 



|a| 



2 



(173) 



for m = 4 




(174) 
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and for m = 2 

D Irnr (^) 



1+f 
2 



if\a\ < 1 
if\a\ = 2 
if \a\ = 2 and d = 3. 



<K 2r {(^) 2 if\a\ = 2andd = 2 



(175) 



Proof: Induction in r, with the statement of the Theorem as the inductive 
hypothesis, r = 1 is trivial because Q^i = and because the statement 
holds for I m i(0). Let r > 2 and the statement hold for all r' < r. The 
inductive hy; 
implies that 



inductive hypothesis applies to both factors Imjr k in Qmr- For = 4, it 



f N) (t) 

1 m k r k \°) 



\a\ 



m k TkH 



K f ~\ a \j{\ a \ + ^-2) 



(176) 



for all t > 0. Recall Eq. ( |i~59|) and Eq. (|l6(]), and 



D a K m Q^ r \t) = A m D a Q^(t) < D a Q^ r \t) 



< dK^)ii\\X ati (t)\ . 



(177) 



Let m > 6. By Eq. ( |15SD and mi + rri2 = m-\-2i, the i-dependent factors 
in X a> i(t) are 



e t(HL-2+|ai|) e i(^-2+|a 2 |) e t(2-i+| Q0 |) = e *(f -2+|q|) 

Since ?r — 2 + |a|>^ — 2 > 1, integrating the RGE gives 



(178) 



zr/W(t)| o < |otW( )| q + 1 I 



1 7TJT" I 



1 



m — 4 



+ - I ds 

o 2 

o 

2- 



where 



iW = AB J dRi i!(32J 1 ) i - 1 ^ -^K miri K ri 



(179) 



(180) 



For m > 6, m - 4 > f . Moreover, ^ a < 3 |q| , so Eq. ( p72| ) and Eq. 

(|l73; ) follow. 
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Let m = 4. One of mi and m,2 must be at least six because the ladder 



part is left out in Imr ■ 
is an extra small factor et = eoe - *, so 



Thus i = i(mi + m 2 - 4) > 3. By Eq. ( |l60|) , there 



(*) 



l + t i_| a] (1) 
o 2 



(181) 



which proves Eq. (|7g). Eq. dl74|) follows by integration. 

Let m = 2. The case mi = m-2 = 2 is excluded since this is the skeleton 
RG. Thus i = \ (mi + m,2 — 2) > 3, and Eq. ( 160 ) implies 



< 



l + t 



'd,2 



^2/ 



(182) 



Thus Eq. ( |172| ) holds, and Eq. (p.75|) follows by integration. 



Proof: [of Theorem 



Eq. O, < ^2,r6 e 

F (2V) 



Let \a\ < 1. By Eq. (|175| ), 

as t 



-t l+t 

2 



< K 2 , r . By 
oo. Thus the limit 



t — > oo of (*) ex ists and is a C 1 function of (pi,P2)- All constants are 
uniform in 0. The second derivative of Q^) is 0(t) in d 



2 and O(l) in 

d = 3. Since Q^r = for i > log/3eo, the integral for I^y over t runs only 
up to log/3eo, which gives the stated dependence on log/3eo. ■ 



r ~(7V) 



Remark 12 iVote that the bounds for m = 2 and \a\ = 2 in Theorem ^ 
do not imply convergence. This is the source of the logarithmic behaviour 



discussed in Section 2.6 and in M 



The following theorem implies Theorem [|. Here the bounds have an 
explicit /^-dependence. One could avoid this /3-dependence by including 
polynomials in t, but this would also increase the combinatorial coefficients 
by factorials. 

Theorem 7 Let a be a multiindex with \a\ < 2. For all m,r, let L mr (0) be 
C 2 in (pi, . . . , p m ), and assume that there are Kml > such that \D a I mr (0)\ < 
Kml, with Kml = for m > 2r + 2, and K^\ = 5 m ^v, where v > 0. Let 
I mr (t) be the Green functions generated by the skeleton RGE, as given in 
Definition^, with initial values I mr (0). Let K m \ = and for r > 2, let 
K mr be given by Eq. (|171|) . Then for m > 4 



q rn I I 

D a Q mr (t) o <2K mr (log Pe ) r - 2 e t 2 11 (183) 



49 



and 



For m = 2, 



D a I mr {t) <K mr (log f3e ) r - 1 e t 2 H 



D a Q 2r (t) q < 2K 2r (log/3e 



,r-2+5 d ,2 6 2-|a| 



and 



D a I 2r (t) Q <K 2r (log/3e ) 



, -2 J 1 */ |a| < 1 

(log/3e ) 1+(5d - 2 if\a\=2. 



(184) 
(185) 

(186) 



Proof: The proof is by induction in r, with the statement of the theorem 
as the inductive hypothesis. It is similar to the proof of Theorem pi with 
only a few changes. Note that mi = 2 and m 2 = 2 never appear on the 
right hand side of the RGE because of the skeleton truncation in Definition 
0. For m > 4, use Eq. (11591); this gives 



D a Q mr (t) <9u mr (log/3e ) 



-2 2-f -[«] 



(187) 



For m > 6, the scale integral is as in the proof of Theorem ||. For m = 4, the 
scale integral is now J * ds = t < log /3eo- This produces the powers of log @eq 

upon iteration. For m > 2, use Eq. (|160| ); this gives K ^^^(^{1 +t)) Sd ' 2 

in Eq. (j!87|). The theorem now follows by integration 



instead of e t 2 a 
over t, recalling that the upper integration limit is at most log/3eo- H 

Proof: [ of Theorem || Convergence of the selfenergy follows for \a\ < 1 
as in the proof of Theorem ||. For m = 4, the function is bounded uniformly 
in t. For \a\ = 2, convergence at (3 > holds because J^Zir^) = for 
t>log/3e . ■ 



7 Conclusion 

The determinant bound for the continuous Wick-ordered RGE removes a 
factorial in the recursion for the Green functions. If the model has a propa- 
gator with a pointlike singularity, power counting bounds that include this 
combinatorial improvement can be proven rather easily. The improvement 
may lead to convergence, but I have not proven this here. If it does, then 
Theorem [l] implies analyticity of the Green functions where all relevant 
and marginal couplings are left out in a region independent of the energy 
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scale, and Theorem || implies analyticity of the full Green functions for 
\\\(3 d < const in many-fermion systems for all d > 1. Natural models 
to which Theorem [l] applies are the Gross-Neveu model in two dimensions 
and the many-fermion system in one dimension. In both cases, I have only 
given bounds where the marginal and relevant terms where left out, to give 
a simple application of the determinant bound derived in Section ||. In both 
cases, analyses of the full models, including the coupling flows, have been 
done previously (0, [Lj and |§ |17|, |18|). 

Many-fermion models in d > 2 are the most realistic physical systems 
where the interaction is regular enough for the analysis done here to apply 
directly (i.e., without the introduction of boson fields). I have defined a 
criterion for Fermi liquid behaviour for these models. A proof that such 
behaviour occurs requires more detailed bounds and a combination of the 
method with the sector method of p0| . This may be feasible by an extension 
of the analysis done here. The Jellium dispersion relation E(k) = k 2 /2 — fi 
is the case where the proofs are easiest because £|s is constant. The proof 
that such a system is a Fermi liquid is possible by a combination of the 
techniques of |2(| and [^IJ; is may also be within the reach of the continu- 
ous RGE method developed here. Perturbative bounds that implement the 
overlapping loop method of [21] in a simple way were given in Section [6|. 

The verification of regularity property (2) for nonspherical Fermi surfaces 
is not as simple, but the perturbative analysis done here can be extended 
rather easily to include the double overlaps used in [23], because the graph 
classification of [23] arises in a natural way when the integral equation for 
the effective action is iterated. Multiple overlaps can also be exploited in the 
RGE; this may be necessary for the many-fermion systems in d > 3. The 
split of the four-point function into the ladder and non-ladder part done 
in Section 3.2 singles out the only singular contributions to the four-point 
function and the least regular contributions to the self-energy in a simple 
way. The treatment of these ladder contributions was done in |22j| , where 
bounds uniform in the temperature were shown for the second derivative in 
perturbation theory. 

The results in Section 6^2 provide another proof that only the ladder 
flow corresponding to Figure ^ leads to singularities and hence instabilities. 
It should be noted that this statement depends on the assumptions stated 
in Section ^L^, in particular on the choice of eo, in the following specific way. 
The curvature of the Fermi surface sets a natural scale which appears via 
the constants in the volume bounds. Above this scale, the geometry of the 
Fermi surface provides no justification of restricting to the ladder flow. This 
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is of some relevance in the Hubbard model, where two scale regimes arise 
in a natural way: if et > jl, where jl = /i — td, t the hopping parameter, 
is defined such that jl = is half-filling, then the curvature of the Fermi 
surface is effectively so small that one can replace the Fermi surface by a 
square. More technically speaking, the constant Qy, which depends on the 
curvature of the Fermi surface, diverges for jl — > 0. Thus, for small jl, et has 
to be very small for Qyet < 1 to hold. If Qyet > L the improved volume 
estimate does not lead to a gain over ordinary power counting. Only below 
this energy scale, the curvature effects of the volume improvement bound, 
Lemma |(| imply domination of the ladder part of the four-point function. 
To get to scale jl, one has to calculate the effective action. Needless to say, 
the effective four-point interaction at scale jl may look very different from 
the original interaction, and RPA calculations suggest that the antiferro- 
magnetic correlations produced by the almost square Fermi surface lead to 
an attractive nearest-neighbour-interaction [37]. However, one should keep 
in mind that for the reasons just mentioned, it is an ad hoc approximation 
to keep only the RPA part of the four-point function above scale jl, and that 
a correct treatment must either give a different justification of the ladder 
approximation, or replace it by a better controlled approximation. 

Above, I have not discussed how one gets from the skeleton Green func- 
tions to the exact Green functions. The key to this is to take a Wick ordering 
covariance which already contains part of the self-energy. That is, the Gaus- 
sian measure changes in a nontrivial way with t. This can be done such that 
all two-legged insertions appear with the proper renormalization subtrac- 
tions, and it gives a simple procedure for a rigorous skeleton expansion. 
Moreover, it makes clear why it is so important to establish regularity of the 
self-energy: once the self-energy appears in the propagator, its regularity 
is needed to show Lemma ^| (or its sector analogue), on which in turn, all 
other bounds depend. The condition ko > d also enters in this lemma and 
seems indispensable from the point of view of the method developed here. 
Details about this modified Wick ordering technique will appear later. 

There is a basic duality in the technique applied to these fermionic mod- 
els. The determinant bound has to be done in position space, but the regu- 
larity bounds use geometric details that are most easily seen in momentum 
space. The continuous RGE shows very nicely that those terms that require 
the very detailed regularity analysis are very simple from the combinatorial 
point of view and vice versa. 
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A Fourier transformation 

Recall that ifi is defined on the doubled time direction T2 and obeys the 
antiperiodicity with respect to translations by (3, Eq. (|7|), and that n T was 
chosen to be even. T2 is the set T2 = e r Z/2/3Z = T U (T + 0) where 

ft fi 

T = {t G T 2 : r = e T k, k G {-- ^, . . . , -f - 1}}. (188) 

The dual to T 2 is = |Z/2n T Z = {uj = ^k : k G {-n T , ...,n T - 1}}. The 
Fourier transform on T 2 is 

/» = e T £ e- w /(r), /(r) = ^ E e " T /»- ( 189 ) 
If /(r — /3) = — /(r), then /(w) = if ^ is even. In that case, with 

^ t ) = t? E e " r /H ( 19 °) 

with M nr given by Eq. ([H|). The orthogonality relations are 

f dT e iK±.»)T = ^ mn5 1 ^ e w = l (Jr0 _ ^ (191) 

Lemma 9 Let (3 > 0, Eq > 0, and ra T > 2(3(Eq + £? max ), w;/iere -E max is 
defined in Eq. (|20|), and /ei a) 5e defined as in Eq. (|l~4|). T/ien /or all k £ 6 
and a// a; G M nrJ | icD — -E(k)| < E'o implies \u>\ < ^Eq and \E(k)\ < 2Eq, 
and 

i E K(\iQ-E(k)\<E ) <E l(\E(k)\<2E ) . (192) 

ojSM„ t 

Proof: By Eq. (|T4|), Im a) = — (1 — cos(u>e T )) and Re Q = j- s'm(uje T ). 
The condition | icD — £?(k)| < Eq implies |Re u)\ < Eq and |Im u)+I?(k)| < Eq. 
Thus |Im a) I < Eq + E mSLX . Since 1 — cosx > ^x 2 , (Eq + -E max )e r > 
e T |Im Q\ = l-cos(we T ) > ^(ue T ) 2 , so \ue T \ < Ti((3(EQ+E rQax )(2n T )^ 1 ) 1 / 2 < 
5. Since ^ is decreasing on [0, £], £ > |Re Q\ > Ml^^l > ~M- So 
M < Since 1 — cosx < \x 2 , 

7T 2 1 

|Im u>\ < ^{uJS T f < —(3E 2 q— < Eq. (193) 
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Thus |Im u; + S(k)| < Eq implies |-E(k)| < 2Eq. In terms of indicator func- 
tions, this means that l(\iQ-E(k)\<E Q ) < 1(\uj\<^E ) l(\E(k)\ < 2E ) . 
The summation over M„ T gives 

7i. r 

\ l(k|<|So) =\ £ l(|2n + l|<^) <E . (194) 

Note that the last inequality holds in particular if £^ < 1, because then 
the sum is empty, hence the result zero, because 2n + 1 is always odd. ■ 



B Wick ordering 

Recall the conventions fixed at the beginning of Section |3|. 

Definition 4 Let Wrfa.VO = e^r-^Ctfr ^ an & \ et be ^ e Qrass- 
mann algebra generated by (?/>(x)) xe r- Wick ordering is the C-linear map 
ttc '■ «4r — > A' r that takes the following values on the monomials: ric(l) = 1, 
and for n > 1 and X\ , . . . , X n 6 T ; 



n c wx 1 )...i;(x n )) 



n 



Wr(r7,^) 



(195) 



ri=0 



Theorem 8 Wrfa,^) = {e M)r ) ■ Let oti, ... , a n be Grassmann vari- 
ables, and let a{X) = Y^k=i a k Sr(X,Xk). Then 



n c {^{x l )...^{x n )) 



n^)^(^) 



(196) 



Q = 



Proof: By Taylor expansion in 77, and by definition of CIq, 

n £) 

n>o Xi,...,x n er 

dX x ... dX, 



Wr(r/,V) 



yl y 

^— ' n! ^ 

n>0 

V- 



,fc = l 



a=0 



n>0 



,fc=l 



<Sa(*ft) 



Wr(a,V) 



a=0 
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n>0 n ' ^ Vfc=l / 

~ } n c (fa, VOr") = (e ( ^ )r ) . (197) 



n>0 



Eq. ( |196| ) follows directly from the definition of Qc- B 

The next theorem contains the alternative formula Eq. (|48|) for the Wick 
ordered monomials used in Section 0. 



Theorem 9 Let A Ct be defined by Eq. (|40j). Then 

n Ct (?P(X 1 )...iP(X n ))=e- Ac t rP(X 1 )...^(X n ). (198) 
In particular, if Ct depends differ entiably on t, then 

^lcM{Xl) . . . HXn)) = -^pfioWXl) • • • Ip(Xn)) (199) 

Proof: For any formal power series f(z) = E fkZ k '', 

f(A Ct ) e^O* = f Qfa ,Ci r,)A e^, (200) 

so e~ Ac t = e -^Ctv)r+(v,4>)r ^ from which Eq- foH ows by defini- 

tion of £7c t , because derivatives with respect to 77 commute with Ac t . If Cj 



depends differentiably on t, Eq. ( |199| ) follows by taking a derivative of Eq. 
^) because Ac t commutes with 9 ^ c ' . ■ 



C Wick reordering 

m(ri) m(r , 2) 

By Eq. © and Eq. ©, Q r (i,^) = | E EE ^n.mi.ra.ma with 

r iSiifa>i >jii=1 ni2=l 

T 'l+ r '2= r 

= /" cyr G mm (t | ) G m2r2 (t | x( 2 ) ) p(x, ip) (2oi) 
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where = (X u ...,X mi ), X {2) = (X mi+1 , . . . , X mi+m2 ), X = (X 1 , . . . , 
X mi+m2 ), P{X,ij>) = J dXdYC t (X,Y) u x>Y (il>), and 



mi 



m2 



. fci=i 



vfc2 = l 



(202) 

Only even mi and m2 contribute to the sum for Q r (t, ip) because Q r (t, Y>) is 
an element of the even subalgebra for all t. Let 

mi ni2 

V (1) = J>* $r(X,X k ), v ^ =Y^ rimi+k 6 r (X,X mi+k ), (203) 



k=l 



k=l 



and r] = r/ 1 ) + rf 2 \ With this, 



mi d 



(-1)' 



mi 



11 a% <ty(X) 



r)=0 
??=0 

(204) 



?7=0 



since (— l) mi = 1. Similarly, 



1712 



\k=i 



TT 8 ( l) r? ( 2 )(y) e ('? (2) ^)r-2('? (2) ^«'? (2) )r 



?j=0 

(205) 

Since i/ 1 ) is independent of r/ mi +i, . . . , r] mi+m2 , the derivatives with respect 
to r]^ can be commuted through so that 

Z(r],ip) 



dr)i... dr] mi +m 2 



(206) 



J T]=0 



with 



Z(r,,i,) = faW C t r/ 2 )) e^r-K^.D^Wv-K^),^^)^ (2Q7) 
Wick ordering of 4 r! ''^ v , the antisymmetry of D t , and Ct = — Z>t, give 



%^) = (-ie^'^ (2) ^ ft (e<**>r). 



(208) 
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Lemma 10 Let A and B be elements of the Grassmann algebra generated by 
(r](X))x£r , that is, A = ^/ cr a-i V 1 an d B = Ylicr f] > where aj, bj G C, 
and let F(x) be the formal power series F{x) = ^ r >o fr xT - Then 



A A B(r,) F{M)) 



(209) 



where ^ is the derivative with respect to r\, acting to the left. 

Proof: The proof is an easy exercise in Grassmann algebra and is left to 
the reader. ■ 



Let mi + m2 = jU. By Lemma |1C 



_drji ...drjfj. 
By definition of t/ 1 ' and r/ 2 ) 



Z(r},ip) 



Tj=0 



d L d L 



drftV *drp) 



E E 

fcl=l fc2=mi+l 



(210) 



r?=0 



dVk 2 



(211) 



Every derivative can act only once on rji ... rja without giving zero because 
the ?7's are all different. So 



r/i . . . ry M exp 



( 



0, 



LCM1XM2 



■% n 

(fci,fe 2 )6i 



<>L D t (X kl ,X k2 )^)(212) 



drjki 



where Mi = {1, . . . , mi} and M2 = {mi+1, . . . , //}. L = contributes to the 
sum, but gives the t- independent result rji ... 77^, so the i-derivative removes 
this term in Eq. ( 208 ). Let iri(ki,k2) = k± and ^2(^1,^2) = &2- For a term 
given by L ^ to be nonzero, tti\l and 7T2|i must be injective, because 
otherwise a derivative would act twice. Thus the sum can be restricted to 
the set 



C = {L C Mi x M 2 : L / 0, and vr fc | L injective for fe = 1, 2} (213) 
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If Le£, tti(L) =vr 2 (L). Thus 



min{mi,m,2} 

C- |J |J £(B!,B 2 ) (214) 



1=1 SiCMi 
|Sll = |S 2 | = 



with 

B 2 ) = {L G £ : 7n(L) = Bi and tt 2 (L) = £ 2 }. (215) 

Let 

Bi = . . . , 6j}, 1 < 6i < . . . bi < mi, 
B 2 = {fei+i, • • • , 6 2 i}, mi + 1 < 6j + i < . . . b 2i < mi + m 2 , (216) 

then for any L € £(Bi, B 2 ), there is a unique permutation 7r G such that 

l^^A^) : ( 217 ) 

Thus the sum over L splits into a sum over i > 1, a sum over sequences 
6 = (&i, . . . ,6 2 j) with 

1 < bi < b 2 < • • • < h < mi < b i+ i < . . . < b 2i < m x + m 2 , (218) 

and a sum over permutations it £ Si. Therefore 

min{mi,m2} i 

= E EE n^'^w) H{ ?M (219) 

1 = 1 7Tg5i b fc = l 



with 



H(i, 6, tt) = m . . . v , n ( ^- J • (220) 

The derivatives give (since mi and m 2 are even and hence (— l) mfc = 1) 



2 



i(i±l)_ E bk 



H(i,b,7r) = (-1) 2 *=i e(7r) f] % (221) 



fc — 1 
k(2{i>l,...,6 2i } 
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The remaining derivatives acting in Eq. (210) come from the Wick ordered 
exponential in Eq. ( |208| ). They now act on H. By Lemma 10, they give 

8 



k=l \ / 



feg{6l,...,i>2i} 



ri=0 



n 

/ 



fc=i 

,feSf{6 1 ,...,6 2i } 



9% 



> 07,^0 



mi+m2 



1 [ i>(x k 



ri=0 



(222) 



fe=i 
k?{b 1} ...,b 2l } 



The final step is to rename the integration variables to rewrite the Wick 
ordered product in the form in which it appears in Eq. (|52]). Before this is 
done, it is necessary to permute the arguments of the G mr (t \ X^) such 
that Xi, . . . jXfy appear as the first i entries. Since b\ < ... < bi, one 
can first permute {X\, ■ ■ ■ , X^) — > (X^ , X\, . . . , This takes b\ — 1 

transpositions and hence gives a factor (— l) 6l_1 . The next permutation 

(X^jXi, . . . , X^-i, Xb 1 +i, ■ ■ ■ Xt 2 ) — > (Xfej , Xb 2 , Xi , . . . ,Xb 2 -i) (223) 

gives a factor (—l) b2 ~ 2 because X^ has already been moved, etc. Thus 

E ( b k~ k ) 

G miri {t | Xi, . . . ,X mi ) = (— l) fc=1 G miri (t \ Xb 1 , . . . ,Xb z ,X) (224) 

withX = (X P1 ,...,X P ) where {1, . . . ,mi}\{6i, . . . , &*} = {pi, . . . , p mi -i\ 
and pi < . . . < Pm x -i- Similarly, 

i 

E Oi+fc-fc-mi) 

G m2 r 2 (t | Xj) = (-1)*=! G m2T2 [t | X 6i+1 , . . . ,X b2i ,X). (225) 

Thus the 6-dependent sign factor cancels, and upon renaming of the inte- 
gration variables, Vk = Xb k , Wk = Xb i+k , etc., and with m = mi + mi — 2i, 
y becomes 

min{mi,m2} 

3^ri,mi,r2,m,2 — 



dXn Dt m $(X k ) Y(X) (226) 
i=i ^ \fe=l / 



with 



^QO = (-1)^ [ dVdWj2G miri (t \V,Xx)G m2r2 {t \W,X2) 



TTGSi \ fc = l / 



(227) 



59 



and X_, X\, X%, V_, and W_ given in Proposition ^. The summand does not 
depend on b any more, so the sum over b, with the constraint Eq. fl218| ), 
gives 

E = 

B 1 CMl ,-B 2 C-M 2 
\B 1 \ = \B 2 \=i 

The sum over permutations 7r gives the determinant of T>^ } . Finally, 

G min (t | V 1 ,...,V i ,X 1 ) = (-i)^-^ Gmiri (t \X 1} V 1 ,... i V i ) (229) 

and 

G m2r2 {t | W 1 ,...,W i ,X2) = (-l)^G mar2 (t I W h ...,W 1 ,X 1 ) (230) 

which cancels the sign factor, and thus proves Eq. ([53]). 

The graphical interpretation of the above derivation is that G miri (t \ X\) 
and G m2r2 (t \ Xj) are vertices of a graph, with the set of legs of vertex 1 
given by ??i . . . r] mi and the set of legs of vertex 2 given by rj mi+ i . . . rj mi+rri2 . 
The operator ( , D t & ®f 2) ) generates lines between these two vertices by 
removing factors of rj in the monomial r/i . . . rj mi+ni2 . The number of these 
internal lines is i, and the remaining m = mi + ra2 — 2i factors % correspond 
to external legs of the graph, i > 1 must hold because a derivative with 
respect to t is taken, so the graph is connected. The rearrangement using 
permutations is simply the counting of all those graphs that have the same 
value. There are ("I 1 ) ways of picking i legs from the mi legs of vertex 
number 1 and (™ 2 ) ways of picking i legs from the ni2 legs of vertex number 
2, which gives K mim2 i. And there are i\ ways of pairing these legs to form 
internal lines. One could also have used the antisymmetry of G m2r2 (t) to 
get this factor explicitly, i.e. to permute to get 

i 

i\ l[D t (V k ,W k ). (231) 
k=i 

instead of the determinant. This is the result one would have got for bosons 
since there are no sign factors in that case. The order of the W k in G m2r2 (t \ 
Wi, . . . , Wi, X2) is chosen reversed to cancel the i-dependent sign factor for 
7r = id. The graph for ir = id is the planar graph drawn in Figure |2[ The 
sum over permutations ir ^ id corresponds to the sum over all graphs with 
the fixed two vertices and i internal lines. In particular, it contains the 



mi 
i 



777-2 

i 



^m\in 2 i 



(228) 
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nonplanar graphs. If one restricts the sum to planar graphs, only the shifts 
7Tfc(j) = j + k mod i, k G {0, . . . , i — 1}, remain, and this is only the case if 
mi = i of ?ri2 = The factor i\ makes the combinatorial difference between 
the exact theory and the 'planarized' theory. 
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